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1. Introduction 

The concept of solitary waves was brought to light in 1834 by the Scottish civil engineer Scott 
Rusell, chasing a single "wave of translation" on a horse along a channel in Edinburgh. Unlike 
ordinary dispersive waves, these non-linear waves do not fade away and they preserve their shape, 
size and speed even they undergo weak perturbations. A major step forward in the conceptual 
understanding of this counterintuitive phenomenon was the discovery of the Korteweg-de Vries 
equation in 1895 which admits solutions exhibiting the features of solitary waves. Moreover, there 
are stronger relatives of solitary waves among the solutions of the KdV equation -the solitons- 
which even survive collisions amongst themshelves. During the celebrated sixties of the past cen- 
tury, a powerful technique was invented, the inverse scattering method, which allowed the solution 
of many "integrable" non-linear partial differential equations in the class of the KdV equations. For 
a recent review on the concept of soliton associated with all these non-linear equations, see, e.g.. 
Reference [1]. 

Our theme in these Lectures, however, is the analysis of the quantum descendants of these 
classical non-linear solitary waves/topological defects in one and two spatial dimensions. In 1961, 
Skyrme [2] discovered that a certain extension of the non-linear sigma model, the so called Skyrme 
model, has both 3D dispersive and solitary waves among their solutions. Because the model at- 
tempted to describe the low-energy hadron phenomenology, and because solitary waves are formed 
from a heavy classical lump of energy, the idea is natural: upon quantization, dispersive waves 
become light mesons -pions- and solitary waves give rise to heavy baryons -protons, neutrons-. 
This bold idea prompted the task of investigating solitary waves in the quantum domain, mainly 
performed in [3], [4], and [5] as far as our Lectures are concerned. The main examples were re- 
ported and the conceptual framework was extraordinarily well clarified in [6]. We insist that there 
are many more authors who contributed to developing this research topic. Some of them took an- 
other approach, and there are many very good reviews in the Proceedings of several Schools, and 
even important old and modern books. Most of the pertinent bibliography is collected in Refer- 
ence [7]. Here, we only cite the papers and reviews with an approach close to ours: the /z/weak 
coupling/semiclassical expansion. 

The plan of these Lectures is to describe a method for computing one-loop fluctuations in 
one-dimensional and two-dimensional topological defects based on the heat kernel/zeta function 
regularization of ultraviolet divergences. The method began in Reference [8] with the calculation of 
the one-loop mass shift to the masses of the paradigmatic sine-Gordon and A 02 kink as a test. In the 
same paper, some of us calculated the one-loop mass corrections to other two- scalar field theoretical 
models with only one field but with insufficient information about the kink fluctuation spectrum to 
apply the conventional Dashen-Hasslacher-Neveu approach. In References [9] and [10], the same 
group generalized the method to provide one-loop mass corrections for two-component topological 
kinks: i.e., the models addressed have two scalar fields and the kinks considered are such that the 
two components of the scalar field are not zero for the kink solution. The kink fluctuations of 
such type of kinks are determined by non-diagonal Schrodinger operators and the only possibility 
for managing the spectral information needed to compute one-loop mass shifts is the heat kernel 
expansion. 

Prior to our work, at the end of the last century, there had been a renaissance in interest about 
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the problem of the quantization of classical lumps. The new impetus came from the subtleties 
arising in the quantization of supersymmetric kinks. Several groups at Stony BrookAVien, [11], 
[12], Minnesota [13], and MIT [14] addressed mixed issues in the problem by studying the im- 
pact made by using different types of boundary conditions - PBC, Dirichlet, Robin-, regularization 
methods -energy cutoff, mode number cutoff, high-derivatives-, and/or performed phase shift anal- 
yses, in connection with possible modifications due to the quantum effects of the central charge 
of the SUSY algebra. There are several valuable reviews, of different character and scope, in the 
literature on these developments: e.g., [15], [16], and [17]. 

The stimulus to our work on the quantization of solitons came, however, from the discovery of 
kinks in theories with two scalar fields living on a infinite line. In this type of model, there are often 
solitary wave solutions such that the two scalar fields are space-dependent for the kink profiles, see 
e.g. [19], [18], and [20]. The knowledge of the spectrum of these two-component kink fluctuations 
is insufficient to profit from any type of Dashen-Hasslacher-Neveu approach. The only possibility 
is to use the spectral zeta function obtained from the heat kernel asymptotic (high-temperature) 
expansion. This framework was precisely chosen in the SUSY kink problem in Reference [21]. 

Kinks are one-dimensional topological defects, but extremely interesting two-dimensional ex- 
tended structures were discovered by Abrikosov in Type II superconductors [22]. The phenomeno- 
logical Ginzburg-Landau theory allowed the existence of magnetic flux lines when applied to this 
type of superconducting materials. Relativistic cousins of Abrikosov strings exist in the Abelian 
Higgs model and were proposed by Nielsen and Olesen in 1973, see [23], as plausible candidates 
as the basic objects in the early string theory approach to hadron physics. More recently Achucarro 
and Vachaspati have discovered even more complex two-dimensional topological defects in the so 
called semilocal Abelian Higgs model, the bosonic sector of electro- weak theory when the weak 
(Weinberg, mixing) angle is | [24]. This model enjoys a symmetry group that is the direct prod- 
uct of two groups: §U(2) 0U(1). The non- Abelian group §U(2) engenders a global symmetry 
whereas the other symmetry generated by the U(l) factor is local (or gauge). Henceforth, Achu- 
carro and Vachaspati christened the topological solitons of this system as semilocal strings. Given 
the important role that these models play in our present understanding of the Standard model, it is 
convenient to address the problem of studying the quantum behavior of these two-dimensional soli- 
tons, sometimes referred to as ANO vortices or semilocal vortices because the vector (gauge) field 
of these solutions is purely vorticial (rotational). This task was successfully accomplished in the 
^ = 2 SUSY Higgs model independently by Vassilevich, [25], and Rebhan, van Nieuwenhuizen, 
and Wimmer, [26]. 

In the bosonic setting, however, without fermions to cancel a good deal of the bosonic fluc- 
tuations, the problem is much more difficult. A good step forward towards this goal was given 
in Reference [27] in which the authors calculate the energy of the fermionic ANO vortex fluctua- 
tions. We profited from our experience with multi-component kinks to compute the one-loop mass 
shifts of ANO vortices with a quantum of magnetic flux in a purely bosonic setting in [28]. To this 
end we used the heat kernel/zeta function regularization method, jumping painfully from one to 
two dimensions. In [29] our calculations were extended to superposed vortices up to four quanta 
of magnetic flux and we attacked the problem of computing one-loop mass shifts to the topolog- 
ical solitons of the generalized Abelian Higgs model in Reference [31]. We summarized all this 
material in the Proceedings of QFEXT05 and QFEXT07 published in [30] and [32]. 
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One might think that this is a very narrow and highly focused subject. This way of thinking is 
not completely true, for two reasons. First, knowledge of quantum field theories with topological 
sectors other than vacuum sectors is not fully settled down, at least at the level of perturbation 
theory around the ground state. Second, study of the quantum fluctuations around topological 
defects is a problem in the kinship of very important physical phenomena, such as the cosmological 
constant problem and the Casimir effect. Vacuum fluctuations (loop graphs) give rise to a non-zero 
constant term in the Lagrangian of the Standard Model. Coupling of this Lagrangian to gravity 
means that the constant term is a cosmological constant induced by the quantum fluctuations of the 
particle fields of an order of magnitude greater than the experimental value of ^10^^. The Casimir 
effect is an even closer physical phenomenon. Vacuum fluctuations also play a central role. Here, 
the idea is to sum the effect of the vacuum fluctuations in the presence of some set up -parallel 
plates, cylinders, spheres- measured with respect to the vacuum. The outcome is the appearance of 
physical forces on the plates emerging from the vacuum. 

Our goal in this report is to present an analysis of the quantum corrections to the mass of 
topological defects developed in different one-dimensional and two-dimensional systems in the set 
of References cited above in a manner as unified as possible. The contribution is divided into two 
separate parts. In the first, we deal with one-dimensional relativistic scalar fields. We explain the 
problem, the method of solution chosen, and the derivation of a compact formula for one-loop 
kink mass shifts in a multi-parametric family of deformed linear 0(A/^)-sigma models. The 
model is a member of this family for only one scalar field: N — \. Our approach is tested in this 
prototypical case, and detailed computations are offered. We also describe the results achieved in 
another member of this family with N = 2 scalar fields having degenerate families of topological 
kinks. The first part ends with an analysis of the kink one-loop mass shifts of the topological kinks 
in the massive non-linear S^-sigma model studied in Reference [33]. Again the model is embedded 
in the family of linear 0(A/^)-sigma models, taking the formal A ^ oo in the case of N = 3 scalar 
fields, a process that ends with a non-linear field theory. Following [34], we compute the mass shift 
using the Cahill-Comtet-Glauber formula, see [5]. 

Part two is devoted to understanding the quantum fluctuations of semilocal strings and Nielsen- 
Olesen vortices. The action of the semilocal Abelian Higgs model is considered when the space- 
time is the (2+ 1) -dimensional Minkowski R^'^. The mix of local and global symmetries, the 
Higgs mechanism in the 't Hooft renormalizable gauge, and the Feynman rules are discussed. 
We then go on to study the very rich moduli space of topological soliton solutions, all of them 
having vorticial vector fields. Numerical solutions of the first-order field equations arising at the 
critical point between Type I and Type II superconductivity are calculated in the case of circular 
symmetry. These planar solitons become strings seen from (3 + 1) dimensions. The next task is an 
analysis of the semilocal self-dual vortex fluctuations and the subsequent vortex Casimir and mass 
renormalization energies. These ultraviolet divergent quantities are regularized via the spectral 
zeta function of the second-order fluctuation operators as in kink cases. Unlike kink cases, the 
pertinent differentials are 6x6- Matrix second-order PD operators. In the background gauge there 
are fluctuations corresponding to the Higgs field, a real Higgs ghost field, a complex massless 
scalar field, and the two polarizations of the massive vector field. There are also fluctuations of 
the Faddeev-Popov ghost field that restore unitarity by compensation of the real Higgs ghost field 
fluctuations. The heat kernel expansion allows us the calculation of the one-loop vortex mass 
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shift in terms of Seeley coefficients and incomplete Euler Gamma functions. Finally, we provide 
numerical results that suggest the breakdown of the classical degeneracy of the semi-local vortices 
in favour of the embedded Nielsen-Olesen vortices. 



2. Scalar field models in a line 



In this first part we report the procedure of computing the one-loop mass shifts for one- 
dimensional topological defects developed in [8], [9], and [10]. Ultraviolet divergences are reg- 
ularized using heat kernel/zeta function methods - comprehensive reviews of these techniques are 
[35], [36], [37]- and the models considered belong to (l+l)-dimensional scalar field theory. We 
shall always deal with topological kinks, but in a particular model ofN^l scalar fields we must 
struggle with the problem of studying fluctuations of a degenerate-in-energy continuous family of 
kinks. 

2.1 Deformed linear 0(A/^)-sigma models 

We shall focus on a multi-parametric family of deformed linear 0(A/^)-sigma models. The 
target (isospin, internal, • • •) space is R^. Let Xa,(^^ I52, • • • ,A/^, denote the coordinates of a point 
in R^. The (multi-component) scalar fields are maps from the (1 + 1) -dimensional Minkowski 
space to the target space: Xa{y^) ^ Maps(R^'^,R^). Here 3;^, /i = 0, 1, denote the coordinates of 
a point in the Minkowski space-time R^'^ We shall use the following conventions for the metric 
and volume element: 

, 2 . . dXa dXa dXa ^Xa dXa dXa 

The action governing the dynamics of the deformed hnear 0(A'^)-sigma model is: 



S[Xi,X2,---,Xn] = \j dy^ 



<a—l 



N 



HaUbKbXaXb V. ^aXaXa 

a<b 

We choose a system of units where the speed of light is c = 1, but we keep h explicit because 
we shall work in the framework of /z-expansion. In this system of units the dimension of h is 
mass X length, [h] = ML, whereas the dimension of the fields and parameters are: 

[Xa] = m'2L2 , [A] = [Xab] = M-'l-^ , [m] = [ma] = 

Defining non-dimensional coordinates, fields, and coupling constants 

y'2 tn A 
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the action reads: 



I CL—1 M I 



2 

1 



V«=i J a<b 



Besides the usual m parameter, which sets the length scale of the system, and the A, which sets 
the strength of the isotropic quartic couplings, there are ^^^^^^ non-isotropic quartic couplings Gab 
and N parameters giving possible quadratic anisotropics ^ The rationale behind the choice of 
this family of models is as follows: the set of parameters Gab = cr^ = O^a^b corresponds to the 
linear 0(A/^)-sigma model. In this case, there aioN—l Goldstone bosons, owing to the spontaneous 
symmetry breaking of the global 0{N) symmetry. Goldstone bosons do not exist in (1 + 1) dimen- 
sions, see [38]. Even if we were to start from the linear 0(A/^)-sigma model, the (1 + 1) -dimensional 
infrared asymptotics of these massless fields would generate anisotropic quartic and/or quadratic 
terms, such that no global 0{N) symmetry (or any of its continuous subgroups) would remain. 
Note that no cubic or linear terms in the fields are allowed because the discrete subgroup of 
0{N) generated by the internal reflections (pa —(pa^a, would be explicitly broken. 

2.2 Vacuum fluctuations 

In the parameter range 

1>G^> -oo, > max 1 ( 1 + 2cj^^) , 1 ( 1 + 2Gaa) ]>0,a^b, 1 + 2cj^^ > 



the constant configurations 



y l+2Gaa 

are stable solutions of the time-independent and static field equations 
dV 



dcpa 



N 



Y,^b(pb+ L Gab^i + (Ja-^ 

a<b 



These non-zero constant solutions are thus the 2N classical minima of the system. In the quantum 
domain only the absolute minima 



are the true vacua. Tunnel effects triggered by bounces turn the relative minima ) > V{(p^^'^ ) 

into false vacua [40]. The ground states are built from at least one pair of minima (more than a pair 
of absolute minima could exist) and the Z2 symmetry is broken spontaneously (at least) to 

^ A warning: despite the notation, the possibihty that some might be negative will not be ruled out. 
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Let us denote space-time coordinates in the form xq = / , xi = x and let us consider small 
fluctuations (pa{t^x) = (j)^^'^ + Tia{t,x) around a true vacuum. The action at the quadratic order is 

s^V-.nu-.n.] = l/-/.^ I [^^ (-^+.;) (2.2) 



where 



Ma 



2( 



1 



:4(1-C7, 



are the particle masses. 

The normal modes of these system with a infinite number of degrees of freedom are determined 
in terms of the eigenfunctions of the differential operator: 



v 












JknX 



Vi 



To avoid problems with the continuous spectrum, we choose a ID "box" of very large but finite 
length / = and we impose periodic boundary conditions: /^(x + /) = f^{x). Kq therefore acts 
on the Hilbert space = ©^=i lI(S>^) and the eigenvalues 



Kof?i{x) = (ol{kn)f:{x) 



(ol{kn) = kl + ixl 



2k 



are obtained from wave numbers labeled by the integers. 

Classically the system is tantamount to a infinite numerable set of uncoupled oscillators with 
frequencies given by the eigenvalues of that become quantum oscillators upon canonical quan- 
tization. The free quantum Hamilton 



(2.3) 



is given in terms of the creation and annihilation operators, [SJ(^„),Sc(^m)] = ^ac^mn- the quantum 
disguise of the Fourier coefficients. Note that H^^^ is proportional to fi. In general, the operator 
coming from the 2 — jth-order fluctuation term in the expansion of the classical action is 
proportional to fiK Therefore, result (2.3) is obtained in the first-order (one-loop) of the /z-expansion 
(loop expansion). 

The vacuum state, annihilated by all the destruction operators, is a coherent state, and is an 
eigen-state of the field operator: 



0jy;O)=O,a^c, 0,|y;O) 



1 



1 + 2CJ,. 



|y;0) 



It is clearly a ground state of the quantum system with energy: 

(0;V|i/(2)|v.o)^|^Tr,a4 
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2.2.1 Spectral zeta function regularization, the T^o-heat equation kernel, and the T^o-heat 
trace 

We usually measure the energy of any state in QFT with respect to the vacuum or ground state, 
or, equivalently, we set H^^^ |^;0) = ^^Tr^2^Q |V;0) as the zero-energy level. Important physical 
phenomena such as the Casimir effect or the cosmological constant problem have taught us that we 
must be cautious about using this calibration. In particular, in our system there are other topological 
sectors and it is convenient not to set the energy of the ground state in the vacuum sector to zero a 
priori in order to allow a comparison with the energy of the ground state in the kink sector. 

The problem is that j^Tr jiKq is a divergent quantity that, one way or another must be regu- 
larized. Several regularization methods have been proposed in the literature: cut-offs in the energy 
or the number of modes, high-derivatives, etcetera, see, e. g., [3], [4], [14], [15], [17]. As in Ref- 
erence [21], however, we shall regularize the vacuum energy using the zeta regularization method. 
Instead of computing the -trace of the square root of the Kq operator, we calculate the spectral 
zeta function -the —^-complex power of the Kq operator: 

^ °° 1 ^ .471^ 

Uois) = Tr^2 V L TAn^Y-^s = ^ E('^^a\^) ' ^ ^ C . (2.4) 

The series in (2.4) are convergent only if Re^ > ^ although conventionally they are analytically 
continued to the whole s — complex plane to find the Epstein zeta functions, all of which are mero- 
morphic functions of s, see [36]. The central idea of the zeta function regularization method is to 
assign to the divergent vacuum energy the finite value 

{0;V\ {H^'^y\V;0)^r-'hfi (^)'c^o(^) 

at a regular point ^ G C of I^Kq {s). /i is a parameter of dimension necessary to keep track of the 
right dimensions. 

The analysis of the associated ^o-heat equation kernel 

^ + Ko^ Kk, (x, jS ) = , Kk, (x, y; 0) = 5 (x - y) 

will help us to unveil the structure of {^) as a meromorphic function. Here j8 = ^ is a non- 
dimensional inverse temperature because the dimension of the Boltzmann constant in our system 
of units is [ks] = ML. In terms of the eigenfunctions and eigenvalues of Kq, one can express the 
heat kernel in the form: 



I 

The notation used for the Jacobi Theta function, see e.g. [39], is: 



e 







^^\^^^Y^^2Ki{n+a){z+b)+^^T) . a,b^Q,\ , , T G C , ImT > 



2 
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The "modular" transformation t = /^j8^— ^ = allows us to write the heat kernel in the 



form: 



^_./(x-30|. 



because the Jacobi Theta function is a modular form of weight \ (alternatively, this equivalence 
could be derived from the Poisson summation formula). There are thus two ways of writing the 
^o-heat trace (related by the modular transformation): 



(0|/^/3).£.-^A.„^ 



N 







(0|/ 



N 



I 



N 



The Mellin transform of the first form gives the sum of Epstein functions: 

N A^2 N 



1 r°° Alt Ak 



However, Mellin's transform of the Poisson inverted version 



An '-' 



I 



n——oo 
s-l/2 



nl 



identifies the spectral zeta function as a series of modified Bessel functions of the second type, see 
[36]. Moreover, formula (2.5) shows that there are poles of I^Ko{s) only at the poles of the Euler 
Gamma function r(^ — 



1 1 3 5 7 2j + l 



because Ki^2-s{'^l^aln) are transcendental entire functions, i.e., holomorphic functions of s in C/oo 
with an essential singularity at ^ = 0°. The behavior of the heat trace at high-temperature is deter- 
mined by the asymptotic formula up to exponentially small terms: 







(0|/ 



47rj3 



which also characterizes the behavior of Tr ji e P^'^ for very large /. 
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2.3 Kink fluctuations 
2.3.1 Topological kinks 

Because of the ground- state structure of the deformed Hnear 0(A/^)-sigma model, there are 
other static solutions that are not homogeneous. These classical lumps are "one-component topo- 
logical kinks" (TKl) We shall refer to this type of kinks in this way because: 1) They are topolog- 
ical. Their profiles connect one ground state with another when plotted from x = —oo to x = oo. 2) 
They have only one component of the iso-vector field different from zero. 

In order to find their profile, one simply looks for solutions of the field equations such that all 
the field components except the one accommodating the absolute minima are zero: (j)^^'^ = , Va 7^ 
c. Under this assumption the classical energy can be written in the Bogomolny form: 




3 

m 

Therefore, the solutions 
of the first-order equations 



d(^c^2(vm-V{(^r')) 



0r(x) = ±0rtanh[^(x-xo)] 



dx 



are absolute minima of the energy that solve the static second-order field equations. TKl kinks are 
thus space-dependent solutions that interpolate between the two ground states, reached by the kink 
in different components of the boundary of the spatial line at infinity. They have finite energy. 



V2XJ-00 I '"^^ ' '"^^ 3yrr^V2A 

and their energy density is spatially distributed. Despite these features these classical lumps are 
stable because of topological reasons: they belong to topological sectors disconnected from the 
vacuum sectors in the configuration space. 



2.3.2 Small kink fluctuations 

Our goal is to study these sectors in the quantum domain. We start from the fact that the quan- 
tum descendants of TKl kinks are the ground states in the topological sectors. Small fluctuations 
around TKl kinks, 

(l)a{t,x) = (l)^^^\x) + ria{t,x) 

are governed by the quadratic action 



dria dr]a 
dt dt 



- riaKria 
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where the second-order fluctuation operator is a diagonal matrix of Posch-Teller Schrodinger oper- 
ators : 

f-^^^i-cic'scch\^x) ... ... ... \ 



K- 



V ... ... ... -j^^ii^-CNc-sech^[^x]J 

The bottom of the wells with respect to the thresholds jU^ are respectively: — |jU^ and —Cac = 
-2(l + (7«,) 

Like the vacuum fluctuations, the normal modes of kink fluctuations are obtained from the 
eigenfunctions of the K operator: Kaafni^) = ^ai^)fni^)' We impose periodic boundary conditions 
f^{x + l) — f^{x) to escape the problems of the continuous spectrum, and K also acts on: = 

The spectrum is slightly different for the operator acting on fluctuations along the TKl kink 
orbit: 

^cc = - ^ + K -2^c- sech^ [y x] 
which is summarized in the next Table. 



Eigenvalues 


Eigenfunctions 


£2(0) =0 


/^(x) = sech2[fx] 


ec'(3) = l^l 


/^(x) = sinh[fx]sech2[^x] 




fc(^x) = e'^(3tanh2[^x] - 1 - f /A:taiih[f x] - -^k^) 



Table 1: Spectrum for the Hessian operator acting on fluctuations along the TKl kink orbit 

There are two bound states, one of zero eigenvalue due to the breaking of spatial translation 
invariance by the kink, and scattering states with thresholds at (0) = • the scattering 

eigenfunctions one reads the phase shifts: 

f^{x) ""^i^ ^^-(^xii^cW) ^ g^(yt) = -2arctan-^^^ 

The periodic boundary conditions select the even ground states and force the momenta to satisfy a 
transcendent equation: 

kl + 5c{k) ^ Inn ^ , kn —n = - • arctan ^ ^ , G Z 

Although the solutions form an infinite discrete set, they can only be identified graphically. The 
physicist's loophole is to work with very large / such that the information about the continuous 
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spectrum is codified in the spectral density: 



1 2 

■ + ■ 



The operators acting in the orthogonal directions to the orbit, a 7^ c, are of the type: 



Kaa^-^^lil- Cac ' SCCh^ [y x] 



Cac^2{\-a^^ 



again of Posch-Teller type, but a with non-null reflection coefficient in general. To describe the 



spectrum we define the parameter A = J ^^^a + i- 



Eigenvalues 


Eigenfunctions 




//w = 


[sech(fx)H+2-Vi[-7 + 2fl,-7,i-7+A;i(H 


htanh(fx))] 


el{k) = k^ + ^il 


/^w = 


[sech(^x)]'*^2^'i[5 - ik + A, \ - ik + A, 1 - ik;\{\ - 


ftanh(fx))] 



Table 2: Spectrum for the Hessian operator acting on orthogonal fluctuations to the TKl kink orbit 

The integers j — Q,\,2, - • ■ ;/[A — \\ label the eigenvalues and eigenfunctions of the integer part 
of the i4 + ^ bound states. There are also scattering eigenfunctions. Both bound states and scatter- 
ing eigenfunctions are Gauss hypergeometric functions [42] times some power of the hyperbolic 
secant. 

From this information, we obtain the reflection and transmission coefficients, as well as the 
phase shifts, the transcendent equation for momenta complying with periodic boundary conditions 
(PBC), and the spectral densities: 



Ta{k) 



Y{^-ik+A)Y{\-ik-A) 



T{\-ik)T{-ik) 
da{k)^d+{k) + d-{k) ■ 

knl + 5a{kn) = iTlll 



Ra{k) 



Taik)r{l - ik)r{ik) 

~ r(i-A)r(i+A) 

lm{Ta{k)±Ra{k)) 
4^''''^''[Rc{Uk)±R,{k) 



2.3.3 Spectral kink zeta function regularization 

Let us assume, temporarily, that all the eigenvalues orthogonal to the TKl kink orbit are posi- 
tive: £a{j) > 0, Va. This means that the TKl kink is isolated in the configuration space and stable. 
We expand the small kink fluctuations in a basis in formed by the eigenfunctions of K: 



[Baikn] 



fl^ix)+Baikn)e'''-^''^->'flix 
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We denote here the number of even bound states in the direction of component a as Nh{a) and 
remark that the zero-mode eigenfunction in the direction of the kink orbit c does not enter this 
formula because zero modes only contribute at higher orders in the loop expansion. 
Promoting the coefficients of the expansion to the quantum operators 

[Bl{kn)Mkm)] = Sac5mn , [Bl{j) ,8^(1)] = 8ac8ji 

one obtains the free Hamiltonian in the kink sector, assuming ortho-normality and completeness of 
the chosen basis: 



^(2) = 



7=1 



1 



1 



The ground state in this sector (no kink fluctuation at all) is also a coherent state 

4('t«)|^;0) =0, Va, V;t„ ; ^a\K;0) =0,a^ c, $c\K;0) = ±(l)f''\anh[^{x-xo)]K;0) . 
The kink ground-state energy is divergent 

{0;K\Hi')\K;0) = ^Tr,.K'^ 
but we shall regularize it by means of the zeta function prescription: let us take the value of 
{0;K\ (H^^^y'\K;0) = 2'-'h^l (^^^ Tvi^R-' = 2'-'h^ ^^^'^ 

at a regular point 5" G C of C,k{s). The problem is that, since the wave numbers kn are determined 
by a transcendent equation, there is no way of writing the Tr as any manageable series. We shall 
rely on the less rigorous (physicists) formula: 



N /Nbia) . oo 
a=l \ j=l ^a\J) 



1 



where in the / = 0° limit all the bound states Nb (a) = A^^^ (a) + Nb_ (a) must be accounted for and 
the "sum" over the continuous spectra must be weighted with the appropriate spectral density. 
As in the vacuum sector, there is the ^T-heat equation kernel 

d 



dp 



+ K]KKix,y;l3) = 



KKix,y;0) = 5ix-y) 



solved for very large / by: 

N 



trKKix,y;p) = £ 



a=l 



Ma) 1 

.7=1 



2% 



dk 



where Na{k) is a normalization factor that depends on the wave number. The heat trace reads: 



Tr^2^-^^= / ^ dxtTKK{x,x;l5) = ^ 

-^-2 a=l 



I 

.7=1 



/>oo 
J 00 



which can be used to determine the kink spectral zeta function via the Mellin transform: 

1 



1 poo 
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2.3.4 High-temperature expansion of the TT-heat kernel and truncated zeta functions 

We have met two levels of information concerning the spectra of small fluctuations. The 
spectrum of vacuum fluctuations is fully known and it is possible to analytically determine the 
zeta function regularization of the sum over eigenfrequencies in terms of (famous) meromorphic 
functions. This possibility fails dealing with fluctuations around stable and isolated TKl kinks, 
but knowledge of the scattering data, bound state eigenvalues, phase shifts and spectral densities, 
allows us to obtain formulas for the spectral zeta function in the form of integrals (rather than 
series). 

The situation is less favourable in the following three cases: 

1. There exist one or several zero modes but no negative eigenvalues in orthogonal directions 
to the TKl orbit: £^3(0) = 0, c. The TKl kink is stable but degenerated in energy with 
a ^-parametric family of ^-component topological kinks (TKk), also stable, k is the number 
of non-translational zero modes. 

2. There exist one or several negative modes but no zero eigenvalues in directions orthogonal 
to the TKl orbit: e«(0) < 0, a^c. The TKl kink is unstable and decay to some stable and 
isolated two-component topological kink (TK2). 

3. There exist one or several zero modes and one or several negative eigenvalues in directions 
orthogonal to the TKl orbit. Isolated or degenerated families of TKk stable kinks arise. 

The outcome is the same in these three situations when multi-component stable kinks arise: ^ is a 
non-diagonal matrix Schrodinger operator and the spectrum is unknown . 

To cope with the problem of multi-component topological kinks we write the ^-heat equation 
in the form 

and look for a solution based on the ^o-heat kernel: 

KK{x,y;l5) = CK{x,y\^)KK^{x,y\^) 

where the density CK{x^y\^) satisfies the infinite temperature condition Q:(x,3;;0) = Ia^xa^ and the 
transfer equation: 

{j^ + ''-^^^-J^^^'^a.i^^y'^^)=tuac^^^^ . (2.6) 

Then, we seek a power series solution Ck^^ (x, j8 ) = L^^o (-^^ of (2.6), which is tantamount 
to the recurrence relations 

nctix,y) + {x-y)^ix,y) = -^{x,y) + £ C/«c(x)cf_i(x,y) + {^^ - ^'jcf_,{x,y) . 

(2.7) 
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In fact, only the densities at coincident points x = 3; on the line are needed. We introduce the 
notation ^^^C^{x) = \imx^y^^^{x^y) to write the recurrence relations between these densities and 
their derivatives in the abbreviated form 



I 



to be solved starting from: (^)Cq^(x) = S^^Sab- The n-th density is determined in terms of the 
{n — 1) -order coefficients, their second derivatives, and the potential matrix elements Uabi^)- 

The solution of these recurrence relations is achieved by using a symbolic program run in 
Mathematica. We list the lower densities up to the third order: 



+ ^(Atfe' - Ma ) (t^i? W + 2t/i(x)) + - ^^c)Uac{x)Ucb{x) + (^1 - ^l)Uab{x)^ 

There is an interesting point about these densities: the diagonal components are the infinite con- 
served charges of some matrix KdV equation, see [41]. Consider the family of differential operators 

K{^) = + diag(Ai2, • • • , Ai^) - f/(x,i8), 
where the family of potentials solve the matrix KdV equation: 

+ 3 [/— + —[/ + 



dp \ dx dx J dx^ 
The j8 -evolution of i^(j8) can be expressed in the Lax pair form 

such that it is iso-spectral. The diagonal densities codify the spectrum of K. Ergo, c^^{x^x) ^ /(j8). 
Integration of the densities gives the asymptotic expansion of the ilT-heat trace: 

Tr^2e-^^ = tr/ dxKK{x,x;fi)=Y,^^"LC[KW--^ , C[K] = \im dxC{x) 

2.3.5 Truncated zeta function 

Finally, we split the Mellin transform 

Q{s) = J-^^" (Tr,2e-^^-no) = f^/^^/^^"' (Tr,2.-/^^-no) +BK{s;b) 
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into meromorphic and entire parts by choosing an upper limit in the j8 -integration to be optimized 
in each particular problem. The number of zero modes, no, must be subtracted. Truncation of the 
number of terms kept in the meromorphic part 

provides us with a practical formula for the spectral zeta function to be used in our computations 
in terms of the incomplete Euler Gamma functions: 7[z; c] . The two free parameters b and Nq are 
correlated: the larger /?, the greater the No must be chosen to achieve good approximations. 

The subtraction of zero modes that we have performed is a very tricky affair. We split the 
improper integral into two parts: 



no 

I[no] = h [no] +h[no] = lim -— ^ 
e^O 1 [S) 



/ dl5 P'-^ e-^l^ + / Jj8j3^-i 

Jo Jb 



We neglect I2 [no] and regularize the divergent integral /i [no] for Re^ < by assigning to it the value 
of the analytic continuation of I^[no] = b^^\s)' ^^^i^ Re^ > 0, to the whole complex ^-plane. 

2.4 One-loop kink mass shift formula 

We shall now perform the renormalizations needed to tame the ultraviolet divergences, part of 
which we have already regularized using the zeta function regularization method. 

2.4.1 Kink Casimir energy: zero-point renormalization 

The first renormalization consists of subtracting the vacuum zero-point energy. By analogy, 
we shall call this quantity the kink Casimir energy. In the Casimir effect, the vacuum fluctua- 
tions are subtracted from the zero-point energy around some geometrical set-up, e.g., two infinite 
impenetrable plates. 

The regularized kink Casimir energy is: 

(2 \ 

The proper kink Casimir energy is the value of this quantity at the physical point s = 

AEi = A£f (.) = ^ - ;ko{-1)) . 

Because 5" = — ^ is a pole of both 1^k{s) and I^Ko{s), with different residua in (1 + 1) -dimensional 
models, the kink Casimir energy is still a divergent quantity. 

2.4.2 Kink energy engendered by mass renormalization 

There are other ultraviolet divergences due to one-loop graphs that are controlled by the 
normal-ordered of the quantum Hamiltonian. The contraction of two fields at the same point in 
Minkowski space is the sum of a normal ordered product, all the creation operators placed to the 
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left of the annihilation operators, and a divergent integral that in the normalization ID "box" be- 
comes a familiar divergent series: 



The quantum Hamiltonian 



dk 



1 



1 



1 



a—\ 



is normal-ordered by applying Wick's theorem: 

1 -exp[-^ £ ^Ata'^]^ \-='^ + f^L^^a-^ 
a=l ^Ya / a=l "ra 



a process that is equivalent, at one-loop order, to adding quadratic counter-terms to the Hamiltonian. 
Regularizing the divergent coefficients of these counter-terms by means of the spectral zeta function 
of^o 



ir(5+i)^^ , 



we obtain at a regular point of I^Kq + 1) the regularized contribution to the kink energy due to the 
mass renormalization counter-terms measured with respect to the mass renormalization vacuum 
energy: 

A£|(.) = 5y.l j dx i{0;K\ : : \K;0) - {0;V\ : : |V;0) j 



-lim . ^ 



r(^) 

The expression in the second line of the formula is derived from the fact that normal-ordered 
products of field operators acting on coherent states select the ordinary product of the field config- 
uration characterizing the state, see [6]. Note also that we have regularized the divergent graphs 
using l^Ko + 1) instead of I^Kq (s). The reason for this is the convenience of comparing the residua 
of AE^{s) and AE§{s) at the same pole: s = 

Finally, the one-loop mass shift formula is found by seeking the physical value of the s param- 
eter: 



AMk= lim {AE^{s) + AE§{s)) 



(2.8) 



2.5 N=l: the A(0'^)2model 

We briefly describe the standard A(0"*)2 model. There is only one real scalar field and no 
deformation parameters. The action is: 
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The very well known vacua and kinks of this model, as well as the differential operators governing 
the small vacuum and kink fluctuations are written below: 



=±1 , ^0^ 



-^+4 ; 0-W = ±tanh(x-xo) , ^-"^+4---^ 





Figure 1: (left) 3D plot of V as a function of both and x. The kink is also shown as a black line crossing 
from one vacuum valley to the other, (right) Graphic of the potential well A — U{x) in the Schrodinger 
operator K. 

2.5.1 Kink spectral heat and zeta functions 

Because the spectrum of K is fully known, both the kink heat function and all the Seeley 
coefficients can be computed exactly ^: 



/ / 



+ V4^(/Erf[yj8] + e^/^ErfpV?]) 



£ c„(/C)j3«-i/2 ; 



coiK) = / , c„{K) 



2"+i(l + 22"-i) 
(2«-l)!! 



n > 1 



where Erf [z] is the error function. The exact kink zeta function and the truncated theta functions, 
with the zero mode subtracted, are respectively: 



1 i ^ / 2 ^ J 1 3 I, 1 \ , r 



(2.9) 



1 



r{s) Jo 



s-l 



-413 No 



^c„(^)/3«-'/'-l 



V47rr{s)\ {4by-'2 t\ (2^-1)!! (4/7)^+^-5 



y[z\c\ are incomplete Euler Gamma functions and 2/^i[^,^,c;z] are Gauss hypergeometric func- 
tions. We follow the notation of [42]. 



^In [8], these coefficients are obtained by integration of the densities solving the recurrence relations. 
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2.5.2 One-loop A^'^-kink mass shift 

The limit ^ ^ — ^ is very delicate: it is a pole of C,k{s). To take this limit properly we look at 
the regularized formulae for ^ = — ^ + e and then allow e to go to zero: 



lim 



2 V27r e^o V 
hm 



^ lim 
hm 



e 



2 

3 
4 



(-i + e)4-z+^ 



• In ; - 3 In - ^''''^ [i , 0, | , - i] + o(e) 



^ lim 
2\/2;r e^o 



- + 31n-^ 
£ m^ 



K 



where 2^/^'^'^'^^ [a, /?, c;z] is the derivative of the hypergeometric function with respect to the second 
argument. The same strategy is used with /\E^{s) to find: 



mm flLi^ 
lim I — ^ 



/In £-0 
3hm 



4-^r(£) 
r(-U^) 



^ lim 
2\/2;r e^o 

lim 



-+ln^-ln4+(v/(l) 
3 



V^(-2)) + «(£) 



+ 31n 



2^' 



m^ 



-2(2+1) 



where i//(z) is the digamma function: i. e., the logarithmic derivative of the Euler Gamma Function. 
Therefore, the divergences at the pole cancel exactly and we are left with the finite answer: 



AMk = ^Jim^^(A£f(.) + A^|(.)) = (^^ 



2^/6 n\fl 



hm= -0.471 llS/zm 



in perfect agreement with the Dashen-Hasslacher-Neveu result in [3]. We remark that the cancela- 
tions above and those implicit in (2.8) set finite renormalizations, such as the large mass subtraction 
scheme, by imposing the condition that tadpoles vanish, see Reference [21]. 

From the sign of the correction we learn a qualitative fact about the global effect of the scalar 
boson fluctuations on the kink. The kink energy density is the square of the derivative with respect 
to X of the kink profile, whereas the kink energy is the area enclosed by this curve. The decrease 
in kink energy due to kink fluctuations is tantamount to a small decrease in the area. Therefore, 
the net effect is equivalent to a force from the right and another from the left exerted by the scalar 
bosons on the kink in such a way that the kink profile shrinks slightly. Exactly as in the ideal 
Casimir effect with two plates of infinite area. 

As a test to prove the quality of the high-temperature approximation, we now write the result 
for No^lOmdb^l. See [8] and [7] for full details: 



AMr 



^ 10 y[;2 _ 1 41 \ 

0.19947 ^ycn(K) ^ ' ^ hm^-OAlUllhm 



The answer departs from the right result in the four-decimal figure; a reasonably acceptable ap- 
proximation. There are two elements, however, in formula (2.10) to play with: namely, the length 
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b of the integration interval in the Mellin transform and the number A^o of terms in the series re- 
maining. Figure 2 shows that these two variables are not independent. For a given value of A^o, 
there is an optimum choice of b after which the approximate result (it comes from an asymptotic 
series) rapidly departs from the exact result. There is no point in enlarging the integration interval. 
Conversely, taking b longer forces us to increase A/q to achieve an acceptable approximation. These 
remarks are particularly important for light-mass (less than one) particles. The reason is that in 
such cases we need b to be large because, if not, too much would be neglected in the entire parts. 

i.Or 

0.5 - 

rO.f IJO 

Figure 2: (left) Plot of the exact heat trace as a function of j3 with the zero mode subtracted, (right) Graphics 
of the approximated formula, keeping A/q = 2, A/q = 3, • • •, A/q = 11 terms. 

2.6 N=2: the BNRT model 

Next, we address a model with two real scalar fields. This field theoretical system can be seen 
as the bosonic sector of an ^ — 1 supersymmetric Wess-Zumino model of two chiral super-fields 
dimensionally reduced to (1 + 1) -dimensions (plus a reality condition in the fields), see [19]. In- 
teractions are derived from the simplest polynomial holomorphic super-potential accommodating 
the second super-field. In fact, the super-potential was discovered independently by Bazeia, Nasci- 
mento, Ribeiro, and Toledo (BNRT) [18] directly in the (1 + 1) -dimensional field theoretical model 
that we shall describe. 

The following choice of parameters in the general linear 0(2)-sigma model 

9 9 CJ 3 (7^-1 , , 

(7f = , CJ2 = 1-- , (^11 = 2 ' ^^^^^r~ ' l + (^i2 = 2cj(cj + l) 

leads to the one-parametric family of actions: 

- |*f + i^fc' + (l-2CT(CT+l))<.i'fc- + ^fc-} . 

Quartic and quadratic anisotropics are induced by the positive real parameter a G M^. 

2.6.1 Vacua and TKl topological kinks 

The four classical minima are degenerated in energy for any o G M^. Henceforth, all the 
classical minima are ground states, and there are "two points" in the vacuum moduli space. 
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1. Vertical vacua: 



'la ^ ^2g' 

such that the particle masses are 



= 2 (^Yq:^(^ + ^12) - (1 - <7f )J = 2CT , i4 = 4(1 - cr^^) = 2CT 
2. Horizontal vacua 

01 ' ^2 =0 , ^ 2'^" 8 

such that the particle masses are 

= 4(1 - al) =4 , Ml = 2 + ^12) - (1 - ^1)) = 

Therefore, there are also two kinds of one-component topological kinks, see Figure 5. 
1. "Vertical" topological kinks: 

0f(x) = O , (/)2^W = ±-^tanh[y|(x-xo)] ■ 
The vacuum and kink fluctuation operators are: 



C0Sh2[y^x] 



2. "Horizontal" topological kinks: 

0f (x) = ± ^tanh(x - xo) , (l>f{x) = 
The vacuum and kink fluctuation operators are: 



\ dx^ / \ dx^ cosh^x / 

2.6.2 Fluctuation spectrum of horizontal TKl kinks 

We shall not discuss the vertical sector here; the main results can be found in [10]. And 
it is more instructive to explain the horizontal sector. In the kink orbit direction, we meet the 
Schrodinger operator of the A^'^-kink. There is no need of re-compute again the effect of the 
fluctuations parallel to the orbit; we simply use the results of the previous model. 

Fluctuations orthogonal to the orbit, however, are governed by the Schrodinger operator shown 
above. We shall focus for a while on the case in which <7 = / G Z+ is a positive integer. The reason 
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for this is that in this discrete set of models transparent (reflection coefficient equal to zero) Posch- 
Teller potentials arise: 

^^^^ 1^2 A^+l) 

dx^ cosh^x 

The previously defined parameter A determining the number of bound states becomes A = / + ^, 

2 

and the threshold of the continuous spectrum is: 4% = 7^. The spectrum of K22 is summarized as 
follows: 

1. The eigenvalues and eigenfunctions of the discrete spectrum are: 
4{j) = {2J-j)j , j = 0,l,2,--,J 

^«'W = ^ > ^^W=n;-(-^ + (y-.)tanh.)-i^ , ;>1 . 

There are / bound states with eigenvalues below the threshold starting from a zero mode. 
The highest eigenvalue bound state sits immediately at the threshold of the continuous spec- 
trum. These eigenfunctions are termed half-bound states and they always accompany zero 
reflection coefficients. In the remarkable array of numbers that follows we have collected the 
bound state and half-bound state eigenvalues up to / = 10. 

J=IJ=2J=3J=4J=5J=6J=1J=8J=9J=10 

1 






























3 


5 


7 


9 


11 


13 


15 


17 


19 


4 


8 


12 


16 


20 


24 


28 


32 


36 




9 


15 


21 


27 


33 


39 


45 


51 






16 


24 


32 


40 


48 


56 


64 








25 


35 


45 


55 


65 


75 










36 


48 


60 


72 


84 












49 


63 


77 


91 














64 


80 


96 
















81 


99 


















100 



2. The scattering states are also known and are listed below together with the pase shifts and 
spectral densities. 



&2{k) = 2arctan ^ p-I^P ' 



- -— + ptanhx e 
dx I 



ikx 



From this information, one derives the exact heat function: 
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7=0 



/ r 

271 J -c 



4i I 



as a sum of error functions. The coefficients of the heat function expansion are also easily calcu- 
lated: 



Tr 



-PK 



Cn{K) 



47r „=o 
(2n-l)!! 



I Cn{KW 



co{K) = I 



2n-l 



2^+1 



(2^-1)!! 



j2n-l _^ ^2n{J) -^2^ 

2n 



n > 1 



They can be expressed in terms of Bernouilli numbers B2n and polynomials B2n{J)' 

These systems can be useful as patterns for other systems where no analytical information 
about the spectrum is available and there is no hope of finding the exact heat function. The coeffi- 
cients of the heat function must be similar for potentials with similar thresholds and areas (number 
of bound states). Thus, it is important to know the behavior of the heat kernel coefficients in these 
"integrable" cases. Considering n as a continuous variable u, the function 



fN{u) 



2U+1 



m 



Hp 



2u-\ 



is plotted in Figure 3 for three values of /. The peaks, as they should be, correspond to integer 
values u^n, and the coefficients pass to zero when = 12 for / = 2, = 20, if / = 3, = 30 for 
/ = 4, etcetera. These are the orders at which the double factorial in the denominator dominates the 
numerator, ensuring the convergence of the series despite the fact that the coefficients grow really 
large before these orders for / > 4. 






Figure 3: Plots of the function fj{u)\ (left) J = 2 (middle) / = 3 (right) / = 4. 



2.6.3 Cancelation of divergences: one-loop mass shifts of horizontal TKl kinks 

The consequences of the previous analyses can be seen in the magnitude of the one-loop mass 
shifts. The spectral zeta function 

-PK_ 



1 



/ r{s- 



fj-\ 



+ 



2j 



-J 



^5 (y2_y2).+r'^'f2'2'^''2'/2-/^ 




r(.) 
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is essentially a sum of Gauss hypergeometric functions. For instance, in the case J = 3 



1 



/ r(5-i) 
2 32^-ir(5) 



+ 



2 11 3 _1 11 3 4 1 \ r(s+i^ 

^■2^i^2'2+*'2' 8^ + ^'^ '^2' 2 +^'2' 5^ 9^ j 



one can envisage the general behavior from the array of bound-state energies. The ingredients are 
the bound-state eigenvalues, the difference in the eigenvalues with respect to the thresholds, and 
the thresholds themselves (in the subtraction of the zero modes). It is amazing how cleanly the 
maths capture these fine physical data! 

The kink Casimir energy in the direction orthogonal to the orbit for generic / is: 

■ 2\ ^ 



AE^ = lim I'-'^h 



hm 

2\/2e 
hm 



m 

2\ e 



At(C42(^)-C^022 (5)) 



2m 



2V2;r e^o \ m 



lim 



2\ e 



r(e) 



'J-1 



V 



J 



Simili modo, the kink energy due to the mass renormalization counter-terms of the second particle 
reads: 



hm 



1 



AEI = ^ lim — lim 



I r(5+i 



hm 



lim 



V4i /2^+ir(, 

( 2ju2 ^^ 7 (7+1) T{£) 

J^' 'r(-Ue) 



dx 



/(/+!) 
cosh^x 



2\/2k e^o \ 

For instance applying these formulae to the / = 3 case, we find: 



(r ^ , -21 ^ ( ir(0,lAO)rl r. 3 li , 8\ 



+ 



+ 2^1 [2'"'2'~5i 



5 3 9 
log4-3--log- 



3hm 



K22 



1 2ll^ 
lim - + log - log - - 2 



The divergent terms, as well as the terms depending on the auxiliary parameter ji cancel each other 
exactly , to give: 



hm 



\ (2Ff^'«'«)[i,0,|,-i]+log^) +,F^'°'<')[i,0,i,-|]+log^ + 3 



(2.11) 

Using this result and the analogous one for 7 = 4, we obtain the exact answer for the 7 = 3, and 
7 = 4 systems: 

AMi^=^ = -(0.47113 + 0.766861)^m , AMj^='^ = -(0.47113 + 1.11725)^m 

(2.12) 
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-4 




Figure 4: Potential wells for / = 2, / = 3, and / = 4 plotted together. 



The higher the threshold and the broader the area, see Figure 4, the more negative is the 
correction, because the attraction of the well on the scattered and bound mesons is more intense. 
The high-temperature asymptotic formula also provides very good approximations to these exact 
results. We find: 

AM^£^ = -(0.47137 + 0.76675);im , AM^£^ = -(0.47137 + 1.11723);im, 

taking A/q = 20 coefficients for / = 3 and A/q = 40 coefficients if / = 4 in the asymptotic formulae 
of AM^^^ and AM-^^. The (optimal) truncations in the series have been suggested by the plots in 
Figure 3. For the same reason, we consider A/q = 10 in the truncated series of AM^/=2. 

Even if a is not an integer the one-loop mass shift is still known. The spectrum is more com- 
plex but the main differences passing through / are the lack of the (half) bound state, sitting at the 
threshold, one more bound state and a different spectral density. Because the vacuum half-bound 
state, (the constant function) is not compensated by a kink half-bound state it must be accounted 
for in the vacuum spectral function. The half-bound state weight is ^ (hence the name) owing to 
the one-dimensional Levinson theorem [43], [44]. The jump in the number of bound states mirac- 
ulously conspires with the change in spectral density to produce a quantum correction that is a 
continuous function of <7. Interested readers can find full details about this problem in [10]. 

2.6.4 TK2 topological kinks 

The zero-mode fluctuation orthogonal to the kink orbit is a sign of the existence of other 
topological kinks degenerated in energy with the TKl kink, which must have the two components 
of the field different from zero, thus being TK2 kinks. Fortunately, a fair knowledge of the main 
features and the structure of the TK2 kinks in this model is available. The reason is that the potential 
energy density can be written in terms of a polynomial "superpotential": 





The energy, arranged a la Bogomolny, 




d(t>a dW\ fd^a dW 
dx d^aj \ dx d(\>a 




XV2 



4m^ 



•{W(0,(+oo))-W(0,(-oo))} 



(2.13) 
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shows that the solutions of the first-order ODE system 



dx 



dW 



^ = (-l)«(20? + CT0|-i) 
dx 



a,j3 =0,1 



(2.14) 



are absolute minima of the energy. 

In [19], the kink orbits solving this system of equations were identified for the first time. We 
found the flow lines of gradW between the top and the bottom of W, the kink orbits, in [20] simply 
by taking the quotient of the first equation by the second in (2. 14) to find: 



2^l + o^l-\ 



+ (-1)^ 



d^2 

2CJ0102 



= 



Integration of this first-order equation is achieved using the integrating factor |02| ^- The kink 
orbits are: 



Note that the integration constant c must be lower than the critical integration function c^. The 
reason for this is easy to understand just by looking at Figure 5. The kink orbit joins one 
horizontal vacuum with the other two vertical vacua. Beyond that value the orbits escape to infinity 
giving infinite energy to the associated static-field theoretical solutions. 






Figure 5: 3D graph of — V(0i,02) as a function of (\)\ and 02 (left) Flow-lines: in the ranges c G {—oo^c^) 
(middle left), c = c^ (middle right), and c G {c^oo) (right) 

The energy of all these TK2 kink orbits is the same: 

The zero-energy fluctuations orthogonal to the TKl kink orbit obeys this fact: there is no cost in 
energy fluctuating from one kink in this family to another. 

For generic a, it is not possible go further than this: i.e., knowledge of the TK2 kink orbits and 
their energy. There are two exceptions, however. If c = —0°, we again find the TKl kink profile: 



1 



^l^{x) = (_i)« tanh(x-xo 



(/)i>^W = 



a = 0,1 
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For c = 0, the profiles of the corresponding TK2 kinks are also known analytically: 



(/.fi^W = (-!)«' ^tanh[2(l-CT)(x-xo)] , ^f^^lx) = i^sech[2(l -a)(x-xo)] 

where ai , 0^2 = 0, 1. The two kink orbits are half-ellipses. 

The generic form of the kink profiles can be inferred from the o — \ case. Equations (2.14) 
are separable in parabolic coordinates if <7 = ^, and we can give the analytic expressions of the 
kink profiles [20]: 



(^^^'^\x\XQ,d\ = (-1 



2 cosh((x — xo)) J 

^ 

y^J^ + cosh((x — xo)) 



d^±x 



1 



There are two integration constants that are the parameters of this TK2 family with a clear physical 
meaning: sets the kink center as in TKl kinks. The parameter d, related to the kink orbit label 
c, can be loosely interpreted as the relative coordinate between two basic kinks. 
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Figure 6: Kink profiles corresponding to: (a) d = 0, (b) d = a/03, (c) d = I and f J - 



30. 



In Figure 6, the kink profiles for four values of d are shown. If J = only one component is 
kink-shaped and non-null. For d > 0, all the profiles are two-component; one of the components 
is kink-shaped and the other bell-shaped up to J = 1. If J > 1 one component has the shape of 
two kinks and the other one tends to show two lumps. This Figure is correlated with Figure 7, 



J. 



Figure 7: Energy density ^^[x;0, J] for (a) d = 0, (b) d = ^05, fcj J = 1 and (d) d = a/30. 
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0, (b)d = \ 


/a5, (c) d = 



where the energy density is plotted for the same values of d. It is possible to prove analytically 
that the energy density as a function of x has only one maximum if J < 1 and two maxima if 
d > I. Thus, for J > 1 this parameter is a relative coordinate between the two peaks, although the 
distance between maxima is a transcendent function of d. The meaning of the peaks is also clear: 
these values of d correspond to orbits close to the critical orbit where two TK2 kinks joining the 
horizontal vacuum with different vertical vacua live. These are the basic kinks. All this is explained 
in detail in Reference [20]. 
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For generic values of a, one must rely on numerical integration methods to find the profiles of 
the kink solutions. We do this by solving the first-order equations by standard numerical methods 
with the "initial" conditions: 

The reasons for this choice are two fold: (1) For any kink solution, 0i(x) always has a zero. 
Translational invariance allows us to set the zero at x = 0; (2) To ensure that we will find a numerical 
kink solution, we fix 02(0) on a kink orbit for a given value of a and arbitrary choices of c. 

The numerical method provides us with a succession of points of the kink solution generated 
by an interpolation polynomial. The plots of the numerical results show that the behavior derived 
analytically for the kink profiles when a = ^ is generic. For any value of a, the kink profiles are 
composed of two kinks. The parameter c giving the orbit is related to the separation between basic 
kinks. In some range of c, the two basic kinks completely melt into a composite kink. The precise 
value of c at which this happens depends on the value of O". O" = ^ is singled out, because in this 
case c = is the value where two kinks fuse into a composite kink or, viceversa, a composite kink 
splits into two basic kinks. 

2.6.5 TK2 kink small fluctuations 



We consider now small fluctuations of TK2 kinks: = ^a(^'^c) + T]^(x), where we denote 
rK2 kink profiles in the form (j)]^^^^"^ 
is a non-diagonal Schrodinger operator: 



the TK2 kink profiles in the form (j)]^^^^"^ (x) = 0^(x;c). The second-order kink fluctuation operator 



Kic) 



-^+4-Un{x;c) Unix;c) 



Un{x\c) = -(240f(x;c)+4(j((J + l)02(^;c)-6) 
Ui2{x\c) = i72i(-^;c) = -8(7((7 + l)0i(x;c)02(-^;c) 
U22{x\c) = -{4a{a + l)^^{x;c) + 6a^^i{x;c) - a{a + 1)) 

In Figure 8 the potential wells (or barriers) in the diagonal components of K are shown for the 
integrable case a = ^ for several values of c. The correlation with the kink profiles is clear. There- 
fore, the non-integrable cases have similar wells because the kink profiles obtained numerically are 
similar to the analytic ones. 

The one-loop mass shifts can be computed using the high-temperature formula. We use the 
first-order equations to express field derivatives as functions of the fields (avoiding problems with 
the discreteness of space required in numerical methods) and write the Seeley densities. Then, we 
plug the numerically generated kink profiles into these densities and perform numerical integration 
over the real line, to finally find the coefficients that enter the formula. 

The results are shown in the following Tables for three values of O": 1.5, 2, and 2.5, and 
represented in Figure 9. 
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c =-30 c =-1 




Figure 8: Diagonal components —U\\{x) (red) and —U22{x) (green) of the cr = ^ potential for c = — 30,c = 
-l,c = 0.1, c = 0.245 and c = 0.249. 
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-5 
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Table 3: Quantum Mass Corrections to the TK2 family in the BNRT model with cr = 1.5,2, and 2.5 
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Figure 9: The one-loop mass correction in the cases cr = 1.5, cr = 2.0, and cr = 2.5 



In the O" = 1 .5 and O" = 2.5 cases, the pattern is similar. The one-loop correction is equal in the 
zone where the two basic kinks are fused into a single kink and the degeneracy is not broken at one- 
loop level. For values of c in the zone where the two basic kinks are split, the one-loop corrections 
become more and more negative. The classical degeneracy is broken by quantum fluctuations that 
press the two basic kinks apart from each other. This conclusion could be reached by arguing 
qualitatively from Figure 8. For small c, the potential wells both in the parallel and orthogonal 
directions to the kink orbit are attractive. After c = 1 the orthogonal wells start to become weakly 



-1.2'l 
I 

V ^. 
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repulsive. For larger c, the parallel wells develop two peaks whereas the orthogonal wells become 
strongly repulsive. This is the explanation of the repulsion between the two basic kinks and the 
reason why the less energetic TK2 kink corresponds to maximum separation of the basic kinks, 
even though the classical energy is independent of the distance between centers. The case (J = 2 is 
special. The classical degeneracy does not disappear because a simple change of variables shows 
that this case is equivalent to two independent A^"^ models, such that the one-loop correction is 
twice the correction of the A^"^ kink. 

2.7 N=3: the massive non-linear S^-sigma model 

The last scalar field theoretical model that we shall discuss is the massive non-linear S^-sigma 
model studied in Reference [33]. This is the formal X -^oo limit of the deformed linear 0(3)-sigma 
model, a limit that is only meaningful if the fields satisfy the constraint 

2 

xKt,x)+xl{t,x)+xKt,^)=R' , ^' = T 

such that the target space is a §^-sphere of radius R. The action becomes 

S[XuX2,X3] - ^ / 1 1 1^ ■ If - ochIit,x) - alxi{t,x) - ahiit,x) 

where 2o"f , 2(72 ' ^3 ^^3- 

Despite appearances this is a highly non-linear system due to the constraint between the fields. 
This statement is made evident by solving X3 in terms of Xi and Xi- The action becomes: 

s^y dtdx [da^d^xi + d,x2d'x2 + ^^'^f_!^/!^f^' - ^1 (^^) - • ^2(t,x) } . 

2 2 

The vacua are the North and South poles: ;t;3 = ±1. The parameter < cr^ = < 1 is the 

mass of the pseudo-Nambu-Goldstone boson, the quantum of the ;^2-field. The mass of the other 
pseudo-Nambu-Goldstone boson is 1. 

Interactions, however, come from the geometry: 

- ;^ + ;^ (^1^ + ^2 ) + ^ (^1^ +X2f + "^' {xi d^xi + xid^xi) {xi d^xi + xid^xi) 

which shows that ^ is a non-dimensional coupling constant. 

2.7.1 TKl topological S^-kinks 

Using spherical coordinates in the target space 

=/^sin0(/,x)cos(p(^,x) , ;^2(^-^) =^sin9(^,x) sin(p(/,x) , ;t;3(/,x) =/?cos0(^,x) , 
the field equations become: 

□0 - ^sin20 {d^(pdi^(p- cos^ (p-a^ sin^ 9) = , d^{sm^ edi^(p)-^G^ sin^ sin29 = . 
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On the half-meridians (p — j, (p — ^ and 9 = 0, <p = ;r the field equations reduce to the sine- 
Gordon equation. Therefore, there are two types of two-component topological kinks joining the 
North and South poles: 

1. Kl kinks 

<Pk, (^) = f ' (^) = Y ; ^K, (x) = 2arctan e±'^(-^--»^o) 

Xf^ix)^0 , X^^ix)^±-^^^^^^ , zf'W = ±^tanh[.(x-xo)] 

2. /r2 kinks 

(pK2 (•^) = , (pK* {x) = 7r ; % (•^) = 2arctan ^^(-^-^o) 

= %oshiu-.„)i • ^f'W-" • A:fW = ±«™h[(.-.„)i 




Figure 10: a) Ki and ^2 ((5"^ = ^) kink orbits, b) Ki (blue) and K2 (red) kink energy densities 

In Figure 10 the topological kink orbits are plotted, together with the corresponding Kink energy 
densities. The classical kink energies are: 

and the K2 kinks are heavier than the ^1 -kinks. 



2.7.2 Kink fluctuations and one-loop mass shifts 

The geodesic deviation from the Ki kink orbit plus the Hessian of the potential for the Ki kink 
(in a parallel frame to the kink orbit) reads, see [34]: 



K 



dx^ cosh^crx 



-^ + 1 



dx^ cosh^crx . 



As might be expected, for parallel fluctuations to the orbit we find the same Schrodinger operator 
as for the sine-Gordon kink. What comes as a surprise is that the orthogonal fluctuations are also 
governed by an operator, which is of the transparent Posch-Teller type of the A/^ = 1 class with 
slightly different threshold. 
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The one-loop Ki mass shift is immediately obtained from the Cahill-Comtet-Glauber formula: 

J^j^f^^y t 1 Pt^ 
= y/lmR^G ^[sinViH — sin V2 - Vi cos Vi V2COS V2] + ^(^) 

Here, Vi = arccos(^) = | and V2 = arccos(y). The arguments are the square root of the value of 
the bound state eigenvalue O"^ = 1 — divided by the threshold of the continuous spectrum. For 
instance, as a function of a the correction reads: 



Ei^{G) = V2mR^G- 



whereas for a = ^ we obtain : 



hmo 
V27t 



3hm 



;arccos(c7) 



—) 
6V3J 



'2' V2 2V271 \3 
The second-order fluctuation operator for the K2 kink (in a parallel frame to the kink orbit) is: 



^ _ f dx^ '^^ cosh^x 

I 



'dx^ 





+ C72 



cosh X 



The bound state eigenvalue —G^ in the orthogonal direction to the orbit is negative, telling us 
that the K2 kinks are unstable. The CCG formula applied to calculate the one-loop K2 mass shift 
captures this fact, because V2 ceases to be an angle and becomes arccoshd": 



E^^{G) = Ek,{g)+AMk,{g) + &{'-^: 



n 



Vi = arccos(O) = ^ ? ^2 = arccos(/a) 



= VlmR^ ^ [sin Vi H — sin V2 - Vi cos Vi V2 cos V2] + ^ ( ^ ' 

V2;r o g R^ 



Therefore, 



hmG 
V27r 



^ + \/2-G^-i^G + alog (^\/2-G^-G^ 



has an imaginary part, as corresponds to the energy of a resonant state in quantum physics. We end 
this first part by computing: 



.1. 



m 



hm 



4 + V7 + A/31og 



/?2 



) 



32 



Quantum topological defects 



Juan Mateos Guilarte 



3. Abelian Higgs models in a plane 

The second part will be devoted to discussing the quantum fluctuations of the topological 
solitons that arise in gauge theoretical models in (2 + 1) -dimensions with spontaneous breaking of 
the gauge symmetry. The prototype of these solitons is the Abrikosov-Nielsen-Olesen vortex, see 
[22], [23]. Because the spatial components of the gauge field form a purely vorticial vector field, the 
planar soliton is made from lumps of magnetic field. By extending these objects to the third spatial 
dimension magnetic flux lines arise like those existing in Type II superconductors, predicted by the 
phenomenological Ginzburg-Landau theory. In the relativistic Abelian Higgs model the Nielsen- 
Olesen vortices form strings that might confine magnetic monopoles. The one-loop mass shift due 
to quantum fluctuations of ANO vortices has been computed at the critical point between Type I 
and Type II superconductors in [28], [29], and [30] using heat kernernel/zeta function regularization 
methods. 

It was discovered in [24] that a generalization of the Higgs model, the so called semi-local 
Abelian Higgs model containing a doublet of complex scalar fields, also admits planar topological 
solitons similar to the ANO vortices. The moduli space of these solitons is, however, richer, en- 
compassing both CP^ -lumps and ANO vortices as limiting cases with very complicated mixtures 
of these two classes in between, see [46]. The one-loop fluctuations of self-dual semi-local topo- 
logical solitons have been analyzed in [31] and [32], with some surprising results. Our aim in this 
part is the description of these results. We remark that the semi-local Abelian Higgs model is the 
bosonic sector of the electro-weak theory when the Weinberg angle is | and the Z, gauge fields 
decouple. 

3.1 The planar semi-local Abelian Higgs model 
3.1.1 Action and field equations 

In the semi-local Abelian Higgs model there is a Higgs doublet 



The action is built from the electromagnetic tensor, and interactions between the gauge and scalar 
fields determined from a minimal coupling principle by means of the covariant derivative and self- 
interactions of the scalar field are induced by a quartic potential energy density: 




and an Abelian gauge field: 




1,2 



Lie\]{\) 




The action governing the dynamics of the system is: 



e 
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There are the following (dimensional) parameters: [v^] = M, the vacuum expectation value of the 
Higgs field, the scalar-vector and scalar-scalar couplings, [e^] = [X] = M~^L~^, and the Higgs 
and vector particle masses: = e^v^, ji^ = Xv^. The ratio between masses K*^ = ^ is a very 
important parameter that determines whether we are in a Type I or Type II superconductivity-like 
regime. Needless to say, we continue working with non-dimensional fields, <I>, A^, and Minkowski 
coordinates x^, jU = 0, 1 , 2. The metric and volume integration are chosen such that: 



^juv = diag(l, — 1, — 1) , dx =dxodxidx2 

g2 g2 ^2 

: g^yx^x^ =xl-x\-xl , d^d^ = gi^yd^d^ ^ 5? ~ ~ 



The field equations are: 



3.1.2 Global and local symmetries: vacuum orbit structure 

The SU(2) global weak iso-spin transformations 

O(x^) ^\x^) = exp(-^0 . d)^{x^) 
as well as the U(l) -gauge local transformations 

O(x^) ^\x^) = ^^■«(^')0(jc^) , A^{x^) ^A^^{x^) =A^{x^) + 

are symmetries of this system suggesting the slightly deceptive name for this model. The vacuum 
orbit has a very subtle structure due to the combined action of these two symmetries. 
The manifold of zero energy configurations, ^^{x^) = AJ^{x^) = 0^, is: 

The Higgs vacuum orbit is the unit sphere in C^: the orbit of the point ^q = ^ ^ ^ , the north 
pole of the sphere, under the global SU(2) action, i.e., the Hopf fibre bundle: 

— ^ — > 

given by the action of a U(l) subgroup on each point of the sphere. 

Use of the Hopf coordinates 61,62 ^ [0,2;r], Y ^ [0? f ] the sphere allows us to write the 
SU(2) action in the form: 

such that the Higgs vacuum orbit is parametrized as follows: 

0^ = cos^isini/A , = sin^isini/A , (j)^ = — cos02CosyA , = sin02CosVA 
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Let us consider the one form (0=^ {(l>id(l)^ + (l)^d(l)J) G Q\E>^), The Hopf index, labeUng the 
homotopy class of the third homotopy group of the 2-sphere n3(S^) = Z, is: 



h=^f (OAd(0=^f (j)id(j)2 Ad(j)^ Ad(j), 



Thus, 



(l)id(l)2 Adij)^ Ad(l)J = -cos^eisin^VAcosi/AjyAA A 
and the Hopf index of the vacuum orbit is: 



h = 



7t^ 



2n- / dOicos^Oi - / d\i/sm^\i/cos\i/ 

Jo Jo 



2 








71^ 


4 



= 1 



3.1.3 Higgs mechanism and Feynman rules 

The choice of a point in the vacuum orbit §^ spontaneously breaks the §U (2) -global symmetry 
and one would expect three Goldstone bosons. One of the three Goldstone bosons, however, will 
undergo the Higgs mechanism. We shift the scalar field away from the vacuum in H{x^) -real 
Higgs-, G{x^) -real Higgs ghost-, and (p{x^) -complex Goldstone- fields: 



\+H{x^) + iG{x^) 
V2(p{x^) 



The choice of the Feynman- 't Hooft R-gauge 

R{A^,G) = d^A^^{xn + G{xn , V. = -^/ d'x{d^A^^{xn + G{xnf 

needs a Faddeev-Popov determinant to restore unitarity which amounts to introducing a complex 
ghost field: 

R{A'^ , G') ~ R{A^ , G) + (□ - 1 - H{x^)) ■8a{x^') 

Det|^ = / [dx*{x^)][dx{x^)]^xv\ij dxh*{x^){n-l-H{x^))x{x^)^ • 



The action becomes: 



d^x 



-A^ [-g^^ (□ + l)]Av + d^X^d^'x - tl 



+ \d^Gd^G-\G^^\d^,Hd^H-'^H'^d^,<^*d^i? 

- y^(^^ + G^) ^Ay,{d^HG-d^GH) +H{A^A^' - x*X) + iA^{(p*d^'(p - 99^9* 
+ A^A'^|<p|2-^(//2 + Gy + l(G2 + //2)A^A^-^|<p|2(|,p|2+//2 + ^2^2//) 



•(3.1) 



The Feynman rules are summarized in the next two Tables. Table 4 gives the propagators in the 
/?-gauge. There is a Higgs field propagating with mass k^, a Higgs ghost of mass 1, and a vector 
particle with mass also 1 . A Faddeev-Popov (anti-commuting) ghost of mass 1 must be included to 
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Particle 



Field Propagator 



Diagram 



Higgs H{x) 

Higgs Ghost G{x) 

Complex Goldstone (p{x) 

Ghost x{^) 



Vector Boson 



ieh 



v{k^ - K*2 + ie) 
ieh 

v(yfc2-i + /e) 
ieh 



v{k^ + ie) 
ieh 

v(yfc2-i + /e) 

—iehg^^ 
v(yfe2-i+/e) 



Table 4: Propagators 



y^rto Weight Vertex Weight Vertex Weight 



Vertex Weight 




-iK^^ 



'he 





he 



-IK' 



•2 



" 2i—9^^ 



■2_L 
he 



-IK 



2v_ 
he 



-IK 



.2 V 



Table 5: Third- and fourth- valent vertices 



compensate the non-physical Higgs ghost. Finally, there is a complex Goldstone boson, as expected 
from this partial Higgs mechanism. The vertices are read from the cubic and quartic terms in action 
(3.1) and shown in Table 5. There are accordingly three-valent and four- valent vertices. The cubic 
terms with derivatives in (3.1) provide weights proportional not only to the coupling constants but 
also dependent on the momenta. 
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3.2 Topological solitons 

The search for time-independent finite energy solutions requires the use of the Weyl gauge: 
Ao(x) = 0; otherwise, time-dependent gauge transformations would spoil the time independence. 
We thus look for solutions of 

\)-d^Ai + diFij^'-{^^Di^-{Di^)^^) , 2)-^5o'^+^AAcI>=|<I>(^^<I>-l) (3.2) 

such that their energy 

E{^,Ai)^ I j2^|l77.^./^.^. + l(Ac^^ (3.3) 

is finite. 

3.2.1 Topology of the configuration space 

The configuration space 

^=[<t>(x) eMaps{R^,C^),Ai(x) e Maps(R^ ,TR^) /E{<t>,Ai) < +00} 

is the set of all the field configurations of finite energy. Consider polar coordinates in the plane: 
r = +-^2' ^ ^ arctan^. The equations 

1) c&^cDlgi =1 , 2) A-c&lgi = {di^-iAi^)y^=0 (3.4) 

are the necessary conditions to be satisfied by finite energy static field configurations at the bound- 
ary of M? at infinity: Si, = {xi,X2/lim^^oo(x^ = r^)} = dM?. Therefore, equation (3.4(1)) 
determines a map 

from si in If<I>2|g, =0 



iie 



are all the single- valued maps complying with (3.4(1)). Acting with the §U(2) matrices parametrized 
by Hopf coordinates we obtain the general solution by changing the base point in §^ at which §^ is 
fibred: 

_ ( ^'^^siny/ e'^^co^^ \ ( ^'^^ \ _ / ^^(^i+^^)sini// \ 
^ ~ I -^"^^2cosi// ^-^^isini// y I j ~ I -^-^(^2-/0)cosi// ) 

Equation (3.4(2)) is consequently solved by: 

showing that the vector field A/(xi,X2) is asymptotically vorticial. The integer number / has corre- 
lated mathematical and physical meanings: 



37 



Quantum topological defects 



Juan Mateos Guilarte 



1. It is the winding number of the map provided by the Higgs field at infinity: 

where §} is the fiber at the north pole of S^. 

2. It is the magnetic flux of the field configuration: 

f f f dQ do \ f^^ 

g=(h (Ai(x)dx^ +A2(x)dx^) =1 (p { —dx^ + —dx^] =1 dO^lnl 

hi ^ ^ / /gi^ J Jo 

Because the first homotopy group of a circle is non-trivial, 111 = Z, the configuration space is 
the union of numerable infinite disconnected sectors distinguished by the integer /: ^ = 

3.2.2 Self-dual semi-local topological solitons 

At the critical point between Type II and Type I superconductivity where the masses of Higgs 
and vector particles are equal, K*^ = 1, it is possible to write the energy, up to a total derivative, as 
a Bogomolny splitting: 

The solutions of the first-order system of partial differential equations (3.6) 

Di^±iD2^ = , Fi2±^(<l>^^-1) = (3.6) 

are absolute minima of the energy and saturate the topological bound proportional to the quantized 
magnetic flux in each topological sector. Because the first-order vortex equations can be obtained 
as a dimensional reduction of the self-duality equations of Euclidean Yang-Mills theory in four 
dimensions, solutions of the PDE system (3.6) -that also solve the second order PDE system (3.2)- 
are usually called self-dual. 

The structure of the moduli space of solutions of (3.6) has been completely unveiled in [46], 
see also [45]. The parameters underlying the 41 dimensional moduli space of topological solitons 
are the coordinates of the / zeroes of ^i, the coordinates of the / — 1 zeroes of and the scale 
and phase of ^2- Full details can also be found in [31]. 

3.2.3 Self-dual semi-local topological solitons with mixed circle-symmetry 

We shall restrict ourselves to study solutions enjoying symmetry with respect to combined 
circle transformations in the M? plane and the internal space C^. This symmetry is materialized by 
means of the ansatz: 

/ f(r)e^^^ \ air) 
^(xi,X2)= r 1 , /,mGZ+,a)GM , A,-(xi,X2) = -/e.-j^x,- . 

Note that by taking some care with the behavior of a at the origin the vector field A/(xi,X2) is 
divergence-free (purely rotational) in the whole plane. The PDE system (3.6) reduces to another 
first-order non-linear ODE system linking /(r), a(r), and h{r)'. 

(3.7) 
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Finite energy solutions, regular at the origin where the (multi) vortex sits, require us to solve (3.7) 
with asymptotic and core behavior: 

lim/(r) = l , lim/z(r)=0 , lima(r) = 1 (3.8) 

/(0)=0 , |MO)| = |/zo|5^,o , a(0)=0 , m<l (3.9) 

We stress that the contribution to the vorticity of the ^2 field m must be always smaller than the 
vorticity / of the ^1 for the topological solution to have finite energy. For later use, we give the 
magnetic field B{r) = ^ ^ and the energy density of these circle configurations: 

Sir) = \{], + m-f{r) - \h{r)\'f + " a{r)f + _ a{r)f . 

3.2.4 Topological solutions of one quantum of magnetic flux 

We now go on to the most elementary solutions that carry a quantum of magnetic flux, or, 
/ = 1 = m + 1 . We are guided by the procedure developed in [47] to solve the non-linear ODE 
system (3.7) with boundary conditions (3.8)-(3.9). First, we consider small values of r and in the 
first-order differential equations we test the power series 

/(r) = fi-r + f2-r^ + fi-P+f4rU--- (3.10) 
a(r) = ai-r+a2-r^ + a3-r^ + a4-r'*-\ — (3.11) 
h{r) = hQ + hi-r + h2-r^ + h-r' + hA-r^ + ■•■ , (3.12) 

where fj and ccj, j — 1 , 2, 3, • • •, are real, whereas hj, j = 0, 1 , 2, • • •, are complex coefficients. The 
coupled first-order ODE's are solved at this limit by (3.10)-(3.1 1)-(3.12) if 

f{r) ^ /, . r + ^(l^op - 1) • + A [(|/,,|2 _ i^^2\ho\' - 1) +4/f] ■r' + ... 



a(r)c^\il-\ho\')-r' 



^ \ho\H\ho\'-l) + lff 



32 



^\ho\H\ho\' - mho\' - 2) + ^ff{5\ho\'-4) 



h{r)c:.ho+'^{\hof-l)-r' + ^[{\ho\'-l){2\ho\'-l)+4f^]-/ + ... . 

We stress that ho G [0, 1] is determined by the behavior of the solution at the origin such that only 
a free parameter, /i , is left. Second, a numerical scheme is implemented by setting a boundary 
condition at a non-singular point of the ODE system, which is obtained from the power series for a 
small value of r (r = 0.001 in our case). This scheme prompts a shooting procedure by varying /i, 
where the correct asymptotic behavior of the solutions is obtained setting a optimal value for /i for 
a given value of ho. Finally, the first-order ODE system is solved for long r by means of a power 



series in ^: 



fir) =Zf'- «(^) = I • Hr) = I /J^' ' 

;=0 ;=0 ;=1 
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with the result that 

f{r) ~ 1 - . + {-2\h' \h' |4) . r-4 + \h' \\-32 + 5\h'\^-^\h' \') ■ r"*^ + 
z o lo 

+ |/ji|2(_ii52+158|/ji|2-^|/zi|4 + ^|/ji|6).^-8 + ... 
I I V II 4 I I 128' ' ^ 

a{r) ~ 1 - \h\^ ■ + \h^\^{-% + \h^\^) ■ + \h^\^{-\92 + 2A\h^\^ -\h^\'^) ■ 

+ \h^ 1^(9216 + 1 120|/i' p - 48|/z' |^ + |^) • r"^ + . . . 

|/i(r)| ~ |/z|i-r-i + -^-r-3 + |/zY(-2+^|/i^P)-r-5 + 
2 8 

+ |/zi|3(-32 + 5|/ii|2- A|/jl|4).^-7 + ___ 
16 

Again, only one free parameter is left. The value of is fixed by demanding the continuity of 
the solution at intermediate distances (r = 15 in our case) obtained by gluing the short-r and long- 
r approximations. In particular, this has the important implication that |/zo|=0 ^ |/z^|=0 
linking the null value of |/zo|, which gives the embedded ANO vortex, with the null value of the 
constant \h^\ setting the behavior of the solution for very long r. Another important remark is that 
the long r behavior of self-dual semi-local defects differs from the long r behavior of self-dual 
ANO vortices that decay exponentially. 

The following Figures show the results obtained with this procedure for several values of 
/zq. Note that /zo = for the ANO vortices, and /zo = 1 for the CP^ -lumps. It may be observed 
in the graphics that the field profiles reach their vacuum values at distances of the order of r = 
15. Consequently, practically identical numerical solutions would be generated by sewing the 
numerical to the asymptotic solution at r greater than 15. 




a) Functions f{r), h{r) and a{r) and b), c) Energy density for /zq = 0.1 




a) Functions f{r), h{r) and a{r) and b), c) Energy density for ho = 0.3 
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a) Functions f{r), h{r) and a{r) and b), c) Energy density for ho = 0.9 
Figure 11: Field Profiles and Energy Densities for Semi-local Topological Defects 

3.2.5 Deformation of the first-order equations 

It is interesting at this point to consider small deformations of the solutions of the (3.6) system 

^{x) = 5(3?) + 5S{x) , Aj (3c) = Vj (3c) + 5a j (3c) 

which are still solutions of the same PDE system. The necessary and sufficient conditions for this 
are tantamount to the linear PDE system: 

-^25^1 + di5a2 + ^(5^55+55^5) = 

(4- - iVi (3c) + i4- + V2(x))dS - i(dai + i5a2)S = . (3.13) 

OXi OX2 

Pure gauge fluctuations are discarded from the solutions of (3.13) by setting the background gauge: 

djdaj{x) + ^{S^5S-5S^S) = 

The tangent space to the moduli space of self-dual topological solitons with a given magnetic 
charge Inl is therefore the kernel of the first-order deformation operator ^: 



l-di 




s\ 


si 


s\ 


Si \ 




/ 5ai (x) \ 


-0*1 






s\ 


-si 


s\ 




5a2{x) 


S\ 


-s] 


-^2 + ^1 


-di-V2 










dS\{x) 




s\ 


di+V2 


-d2 + Vi 










dSlix) 


s\ 


-Si 








-d2 + Vi 


-di-V2 




dS\{x) 


\ Si 


si 








di+V2 


-d2 + vj 




\SSl{x)) 
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Accordingly, A+ = is the following 6x6 matrix partial differential operator: 

I A -2ViS\ 2ViS\ -IViSl 2ViS\ \ 

A -2V2S\ 2V2S\ -2V2SI 2V2S\ 

^ _ -2Vi^f B -2Vkdk S\S\^SISI Slsj-S^S], 

+ ~ 2Vi^} 2V2S\ 2Vkdk B -S\sl^Sls\ SlS^^S^Sj 

-2ViSl -2V2SI S\s\^Sjsl -S\sl^SlS\ C -2Vkdk 

\ 2Vi^i 2V2^| ^}^^_^2^i ^1^1+52^2 2Vkdk C J 



-dkdk+\Si\^ + \S2\\ J,^=l,2 



B 



-dkdt + -{3\Si\^- 



■\S2\^ + 2VkVk-l) 



C = -dtdt + -{\Si\^ + 3\S2\^ + 2VkVk-l) 



,12 



VjSi = djSi + e^%S^„ , M=\,2 , A,B=\,2 , = -e 
One easily checks that A+ has a supersymmetric partner^: A_ — QiQi'^: 



.21 



1, e 



11 



„22 



0. 



/AO 











\ 


OA 













B_ 


-2VA 








2VA 


B- 








s\s\^s\s\ 




c_ 


-2VA 




Vo s\s\-s\s\ 


s\ s\ + s\s\ 


2Vkdk 


C- 


/ 



B. 



-dkdk + ^{\Si\^-\S2\^ + 2VkVk + l) 



v.4 = vi(x)^+y.(x)^ 



C_ = -dkdk + -{-\Si\^ + \S2\^ + 2VkVk + \) . 



The index of the deformation operator - ind & — dimKer^ — dimKer^''^ - is in this case equal to 
the dimension of KerA+ because dimKer^''^ — 0, A_ being definite positive. 

Because A+ and A_ are iso-spectral up to zero modes the index of ^ can be regularized in 
the following form: 

index^= lim [Tr,2e"^^+ - Trrae"'^^- 
Let us use the heat trace expansion: Tij2e~l^^^ — Ttjie'^'^ En=o <^n[A±]j3", 



/-dkdk + l 

-dkdk + l 



A 



V 

















\ 



-dkdk + l 

-dkdk + l 
-dkdk 

-dkdk J 



Because 

^ 7t 



27rp 



co[A±]=6/2, , c,[A±]^ J dxhT{A±{x)-A) 



^This is an example of hidden bosonic supersymmetry [48]. 



42 



Quantum topological defects 



Juan Mateos Guilarte 



we obtain: 

ind^^-^J Jx^tr(A+(3c)-A_(3c)) = ^Jdx^ + |52p - l) = dx^Fn^Al . 

(3.14) 

In the derivation of (3.14) we have used the vortex equation (3.6) and tr means trace in the matrix 
sense. We find that the number of zero modes is twice the magnetic flux (mod tt), in perfect 
agreement with the number of parameters of the self-dual topological solitons. 

3.3 One-loop correction to the masses of semi-local self-dual topological solitons (SSTS) 
3.3.1 SSTS fluctuations 

Let us now consider time-dependent small fluctuations of the self-dual topological solitons: 

^(x) = S{x) + 5S{xo,x) , Aj{x) = Vj{x) + 5aj{xo,x) , x{^) = 5x{x) 

In order to discard pure gauge fluctuations we impose the Weyl/background gauge condition (the 
R gauge in the topological sector of magnetic flux 1): 

Ao(xo,x) = , dj5aj{xo,x) + -(5^(x)55(xo,x) — 5S^ {xo,x)S{x)) = 

The classical energy up to <^(5^) order (one-loop) of the SSTS fluctuations is: 

where 

5(^^(xo,x) = (5ai(xo,x) 5a2{xo,x) 5S\{xo,x) 5Sl{xo,x) 5S\{xo,x) 5Sl{xo,x)) 

is a file vector assembling the fluctuations of the two polarizations of vector particles, the Higgs 
and Higgs ghost fluctuations, and Goldstone fluctuations around the topological soliton solution, 
i.e., the bosonic fluctuations and the second-order operator determining the small fluctuations of 
these extended objects is precisely K = A+. 

We shall impose periodic boundary conditions on the fluctuations 5(^^(xi,X2) = + 
/,X2 + /), a = 1, 2, • • • , 6, / = mL. Therefore, K acts on the Hilbert space = 0^^^ L^(§^ (g)S^). 
Assuming the ortho-normality and completeness of the eigenfunctions of K, K£,n{^) = £^<^«(x), in 
the sub-space orthogonal to its kernel, one finds the quantum Hamiltonian of the one-loop bosonic 
fluctuations: ^ 

H^^^^hm £n{BlBn + -) , (3.15) 

where \B\^Byy\ = bnm are the expansion coefficients promoted to creation and annihilation opera- 
tors. The index theorem argument in the previous Section shows that the number of normalizable 
zero modes is: dimKer^ = 4/ where is the partial differential operator arising in the deforma- 
tion of the first-order equation. Because Qj'^Q), standard supersymmetry strategies allow us to 
conclude that the spectrum of K is formed by 4/ zero modes and positive eigenvalues giving rise to 
Spec^^. 
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The topological soliton ground state is a coherent state annihilated by all destruction operators: 



BnpJS)^Q,\/n 



4>(xo,x)|0;r5) =5(x)|0;r5) 
A-(xo,x)|0;r5)=\^-(x)|0;r5) 



Therefore, the bosonic energy of the ground state in a topological sector of magnetic charge / ^ 
is: 



(0;r5|^PV5;0) = ^Tr^2 7ri 



fim 



(3.16) 



In the Figure 12 we show the diagonal potential wells/barriers in the differential operator K 
for magnetic flux 2k. The vector bosons "essentially" feel -there are also non-diagonal terms- the 
potential A; Higgs bosons and Higgs ghosts feel the potential B, and the Goldstone bosons move 
along the potential C. Leaving apart the non-diagonal exchange interactions, all three types of 
particles move through attractive potentials, exponentially decaying to their vacuum values if the 
background is the ANO vortex, /zq = 0. If the background is a /zq = 0.3 topological soliton, both 
vector and Goldstone bosons move in less attractive potential wells that, moreover, only decay as 
some negative power of r near infinity. For /zq = 0.9, a value giving almost a CP^-lump, the vector 
bosons pass through the background feeling a extremely weak attraction, the Higgs boson and ghost 
note a considerably weaker attraction as compared to the attraction of the previous backgrounds, 
and the Goldstone bosons are repelled by the topological soliton. The decay at infinity is extremely 
slow. The analysis of these physical features will give qualitative support to our results on the 
one-loop corrections to be presented later. 





4 6 S 10 12 14 




Figure 12: Potential wells f or / = 1 VA^r), VB^r), Vc{r)\ a) /zq = 0.0, b) /zq = 0.3, c) /zq = 0.9 

It remains to get rid off the contribution of the Higgs ghosts by considering the ground state 
energy of the (fermionic) Faddeev-Popov ghosts. The FP ghost fluctuations (with PBC) are deter- 
mined by the Schodinger operator: 

such that acts on = L^(S^ Again assuming the ortho-normality and completeness of 

the eigenfunctions of K^, K^Xn{x) = elXn{x), in the sub-space orthogonal to its kernel, one finds 
the quantum Hamiltonian of the one-loop fermionic fluctuations: 

H^^^hm ^ £^(ctQ-l) , (3.17) 

where {C^,C^} = 8nm are the expansion coefficients promoted to fermionic creation and annihila- 
tion operators. There are no fermionic ghosts in the topological soliton ground state. Henceforth, 
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Cn\0; TS) = 0, and the ground state energy of the topological soliton is: 

(0;T5|4')|0;T5) + {0;TS\H^^^\0;TS) = ^ (xr^z/^^ -Tr^2 (/i:^)^) = ^STvi^K^ . 

(3.18) 

3.3.2 Vacuum fluctuations 

It is instructive to specify this analysis for the vacuum fluctuations: 

In this case the second-order operator fluctuations are A^o = A and = —dkdk + l such that the 
spectrum is completely known. In a normalization square of area periodic boundary conditions 
on the fluctuations plus canonical quantization produce the following bosonic and fermionic free 
Hamiltonians: 

a=l ^iGZ n2El^ 
n\ gZ /22^Z 

47r^ 

C0{ni,n2) = C0a{ni,n2) = -^{nl + nl) + l, a = 1,2,3,4 

Therefore, the ground state in the vacuum sector is also a coherent state: 

<t>(^^)|V;0) = 



^ \V;0) 
Ai{t,x)\V;0)=Oi 



such that the ground state energy follows immediately: 

(0;y| (^f +hP) \V;0) = Tr^a/sTj -Tr^2 {l^) ^ = STr^2 (/sTq)^ 

It is clear that the ghosts cancel the contribution to the vacuum energy of the non-physical Higgs 
ghosts and render the theory unitary. 

3.4 Zero point energy and mass renormalizations 

3.4.1 Topological soliton Casimir energy 

Subtracting the ground state vacuum energy from the ground state energy of the topological 
solitons 

AE^s = (0;T5| (^f +Hf ) |T5;0|) - (0; V| (^f +hP) \V;0\) 

Prjfl ( 1 i \ 

= A£ts - A£o = ( STrisTz - STr^J j (3.19) 
one formally measures the semi-local self-dual topological soliton Casimir energy. 
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3.4.2 Mass renormalization counter-terms 

In (2+ 1) -dimensions, the semi-local Abelian Higgs model is super-renormalizable. This 
means that there are a finite number of divergent graphs. We list the one-loop graphs that diverge. 

1. Higgs boson tadpole: 

) = 

= -2i{K^ + 1) • 7(1) - i/c^ • 7(0) + finite part 

2. Higgs boson self-energy: 

+ + s ^ + — < > — + 




= -2/(?c^ + 1) • 7(1) - iK^ • 7(0) + finite part 
3. Higgs ghost self-energy: 

^ Q ^ O ^ O ^ 




= -2/(k-2 + 1) • 7(1) - • 7(0) + finite part 
4. Goldstone boson self-energy: 

+ ('""") + C) + + — r^- 




= -4i{K^ + 1) • 7(1) - HkP' ■ 7(0) finite part 
5. Vector boson self-energy: 




= 2i- [7(1) +7(0)]^'''' + finite part 

Note that the ultraviolet divergences come from integrals of the form: 

d^k i 



7/ 2^ f «^ 



{27ty k^-c^ + ie 

and that, unlike in (1 + 1) -dimensional scalar field theory, normal ordering is not enough. 

46 



Quantum topological defects 



Juan Mateos Guilarte 



Diagram 


Weight 


X 


/[2()c2 + 1)7(1) + ^2^(0)] 


K — 


i[2{K^+\)i{\)+K^m] 


K- 


/[2()c2 + 1)7(1) + jc2/(0)] 


K- 


j[4(/c2 + 1)7(1) + 2jc27(0)] 




-2/[7(l)+7(0)]^''^ 



Table 6: Counter- term vertices 



We add the following counter- terms to cancel these divergences: 

= -;z[/(l)+/(0)].A^(x^)A^(x^) 

We have used a minimal subtraction prescription and the main criteria to set finite renormal- 
izations have been: 1) The divergence due to the tadpole graph is exactly canceled in the self-dual 
limit — \. 2) The global §U(2) symmetry remains unbroken after one-loop renormalizations. 
Interested readers can find a fully detailed description of our renormalization conventions in [32]. 

3.4.3 Mass renormalization counter-term energies 

Therefore, the topological soliton energy due to mass renormalization counter-terms in the 
self-dual limit /c^ = 1 receives the following contribution from the scalar and vector fields: 

A£^r = ^[4/(l)+/(0)] yj2^(l-|5ip-|S2p), A£^r--Sm[7(l)+7(0)]y^2^V,V, . 

We reshuffle the sum of these two quantities into two pieces /^^^ = and 
t\E^^^ = ^ 7(0) (5, Vfc), respectively proportional to 7(1) and 7(0): 

As in the zero point renormalization one must subtract the energy induced by the counter-terms in 
the vacuum from the same quantity for the topological soliton: 

t\E^^ = (7^5) - (^) + A^?^ ijS) - l\E^^^ (y) . (3.20) 

3.5 High-temperature one-loop mass shift formula for self-dual semilocal topological 
solitons 

3.5.1 Spectral zeta-function regularization 

Both t\Ej^ and AE"^^ are divergent quantities that we shall regularize by means of the zeta 
function procedure before being added. We recall that the spectral zeta functions of elliptic op- 
erators A (with positive definite and discrete spectrum) are formally defined as infinite sums of 
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complex powers of their eigenvalues: 

Uis)= £ A- , sec . 

SpecA 

These sums are usually convergent for Re^ > > 0, where is a positive real constant, but are 
susceptible to being analytically continued to the complex 5-plane. In many favorable cases, their 
analytic continuations are meromorphic functions of s and we shall regularize: 1) the ground state 
energy in topological sectors, 2) the ground state energy in the vacuum sector, and 3) the SSTS 
Casimir energy, by assigning to these quantities the values 

A^^^(^) = AEts{s) - AEo{s) , AE^s= 1™ ^Tsi^) 



at a regular point of the spectral zeta functions in the complex ^-plane. The spectral zeta functions 
of Ko and with periodic boundary conditions on the edges of a square of area are given by 
meromorphic Euler Gamma functions: 



n r{s) 2k {s-\)T{s) ' '^0^^ 47r T{s) 

AEjg can be regularized in a similar vein. On a square of area 7(1) and 7(0) become infinite 
sums over discrete momenta: 

1 r 1 _ 1 1 ^ 1 

Therefore, we regularize AE^^ in the same form: 

A^?s(^) = ^(^)'(C-A+i(^)r(^)(5,V,) + C-A(^)r(«)(5,V,)) , AE§s=\imAE^si^) , 
knowing that: 

3.5.2 The heat kernel expansion of elliptic differential operators 

Because the SSTS solutions are not known analytically, it is not possible to compute the spec- 
tral functions of 1^k{s) and {s). One possible loophole is to rely on the high- temperature asymp- 
totic expansion of the heat traces and to build approximations to the spectral zeta functions from 
the Mellin transform of these approximated heat traces. 

This construction is particularly complicated for the second-order fluctuation operator K, The 
Hessian is of the general form: 

K = Ko + Qk{x)dk + V{x) 
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The logical sequence is as follows: the ilT-heat equation kernel is the solution of the K-hQdX equation 
complying with the delta function infinite temperature condition: 

(^^1 + K^Kk{xS^^)-^ , /r^(xJ;0)=r5P)(f-y) . (3.21) 

In our problem I is the six x six unit matrix. The ^-heat trace, or ^-partition function, is the integral 
over the whole -plane of the TT-heat equation kernel at the diagonal of R^(g)R^, whereas the 
spectral TT-zeta function is the Mellin transform of the TT-heat trace: 



Tr^-^^ = tr / d^xKK{x,x;P) , = Tdji P'-'Tre' 

Jr2 i i^s) Jo 



r{s) Jo 

Given the structure of K, 



Kk {x, y;P)^CK {x, f, P)Kk, {x, f, P ) 

is a good ansatz to solve (3.21). Plugging the ansatz in (3.21), one finds the transfer PDE and the 
high-temperature condition that C{x,y;j3) must satisfy: 



(3.22) 



^I + '^^^{dkI-^Qk)~AI + Qkdk + V^CK{x,y;l5) = , C;,(x,j;0)=I 
The solution of (3.22) by means of a high-temperature power series expansion 

oo 

is tantamount to the solving of the recurrence relations: 

[^I+(x^-yy^)(4l-^e)^)]c.(x,3;;ir) = [/\l-Qkdk-V]cn-i{x,f,K) . (3.23) 

It is easy to find the first Seeley density, ci {x^x\K) = —V{x), a result used in the subsection above 
addressing the index theorem. Higher-order Seeley densities are harder to find. It is convenient to 
introduce the notation 

(X) = , [CnU{xrx;K) = (f) 

because the recurrence relations between derivatives of the Seeley densities can be written in the 
compact form 
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-V VV (^^\ f^A P'^'Qf (a,-r+l,a2-t)f^db(^-. , 



susceptible to being solved with the help of a symbolic programm implemented in Mathematica. 

3.5.3 Mellin transform of the heat trace asymptotic expansion 

The spectral zeta functions of of both K and are obtained from the high-temperature ex- 
pansion of the heat traces via the Mellin transform: 

T{s) Ji 



Here the Seeley coefficients [c^^{K)]aa, [c^{K)]aa, and Cn{K^) are obtained through integration over 
the whole -plane of the Seeley densities: 

[ci{K)\aa = J d^x[Cn]aa{x,x;K) , a =1,2,3,4 

[CniK)l, = I d^x[cn]aa{x,x;K) , ^ = 5,6 , [c,(/r^)] = ^ J^x [q] (x ,x . 
Finally, we find the spectral zeta functions as sums of meromorphic and entire parts: 

= t { t ^hsm^ UnmU- 1) } ^ W)''^'^ 

U<s) = £.,(^«)*±^ + J^.,„(.) . (3.24) 
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An important warning: the physical point ^ = — ^ is a regular point -not a pole- of C^(^), 
^^g(^), I^Kq{s), and 1^kg{s). Unlike in kink cases, after zeta function regularization is performed 
only finite renormalizations are left in (2 + 1) -dimensional field theories. 

3.5.4 High-temperature one-loop mass shift formula for self-dual semi-local topological 
solitons 

Neglecting the entire parts and truncating the zeta functions at a finite number of summands 
A/q, the high-temperature one-loop formula giving the semi-classical shift to the masses of semi- 
local self-dual topological solitons is obtained: 



hm 



a^5 ^ 2 



167r0r 

The following remarks are meaningful: 

1. The factor —2/^ is due to the subtraction of the 41 zero modes by exactly the same proce- 
dure as the regularization method used for kinks. 

2. Unlike in theories with only massive particles, the criterion of the exact cancelation of tadpole 
graphs is not completely equivalent to the cancelation of the contributions of the first-order 
diagonal Seeley coefficient. The term in the second row of formula (3.25) is the mismatch 
between these two criteria due to the massless particles in the semi-local Abelian Higgs 
model. 

3. The ANO vortices correspond to the SSTS solitons with S2{x) = 0. Freezing the 552(xo,x) 
fluctuations and dropping away the associated 21 zero modes, the one-loop mass shifts for 
ANO vortices are attained. 

3.6 Mathematica calculations 

In this last Section we shall solve the recurrence relations to find the Seeley densities and their 
associated Seeley coefficients by a mixture of symbolic and numerical Programs implemented in 
Mathematica. We shall focus on circle synmietric SSTS solutions, 

S\{xuX2) = /(r)cos0 5?(xi,X2) =/(r)sin0 

S\(xuX2) = h(r) Sl{xuX2) =0 

(x(r) oc(r) 

Vi{xi,X2) = sine V2(xi,X2) = COS0 , 

r r 

the only ones at our disposal, although /(r), /z(r), and a{r) have been found numerically. 
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3.6.1 Seeley densities for circle symmetric semi-local vortices 

The first step is to use the first-order rotationally symmetric ODE (3.7) to write the field deriva- 
tives in terms of the field profiles themselves: 

rjcos^el M = M [1 -aWsin^el ^ = -^a(r)cose 

dX2 r ^ 6x2 r 

dS^ fir) dS^ h(r) 

a(r) cose sine ^ = -^^a(r)) sine cose ^ = ^a(r)) sine 

oxi r dx\ r 



dx\ 


r ^ 


ds\ 

dx2 


r 


dVi 


cos2e 


dxi 


2 


dV2 
dxi 


= -cos2e 



2/(r)a(r) , /2(r)+/z2(r)-l 



r 2 



dVi ^^Oc(r) sin^e,^9,, ,9,, 

^ = -cos2e^ + ^(/(r) +/i2(r) - 1) 



2/(r)a(r) , /2(r)+/z2(r)-l 



This step is important because numerical calculations on field derivatives cause considerable er- 
rors. The next step is the Mathematica solution of the rotationally symmetric recurrence relations. 
Below we list the three first-order circle symmetric Seeley densities: 



trci(r) = 5-^^-5f{rf-3h{rf 

tr4(r) = ^ [4 a(r)V 27 rVW'- 8 r2a(r) (-1 + 14/(r)2 +Kr)2) + 

+ /(r)2 (^56r^ -64/ + 34/ h{rf^+/ {3I -32h{rf + 1 h{rf 

^'''^'^ ^ T^{-4«W'-4^'«W' (l4 + 35/(r)2-36/.(r)2) + 
+4a(rf (20^9r^ ^32r^ f{rf ^26r^ h{rf^ - 
-2r^a{r) f{rf + f{rf [32 + 33\r^ -15r^h{rf^ -4 (^-\+h{rf^ [-\6-9r^ + h{rf^^ + 
+a{rf [-256 - 144 r^-ni/ + 99 / f{rf - 16 h{rf + 94/ h{rf - 61 / h{rf+ 

+ 2/f{rf [56+\^3/ + \9/h{rfy^ + 
+r4 [- 16 + 15 1 - 29 f[rf + {32 - 135 /^ h{rf + 16 + 23 /^ h{rf + / h{rf+ 
+ f{rf {-2() +199/- 51 /h[rf^+f[rf (^392-321r2 + 2 (-68+ lllr^) /i(r)^ -27r2/z(rf | . 

Here, tr means that we have summed up to the fourth diagonal density. 
trcf(r) = l-^-/(r)2-3Mr)2 

trc2^(r) = ^[Aa{rf-/f{rf + i?a{r) (l +2/(r)2 -^r) 



-^a{rf (l + r^ f{rf -5?h{rf'^ +2? f{rf {-A + 9r^h{r) 
{\-?,h{rf + \9h{r) 
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-4a(r)'^ (20 + 3r2 + 2r2 (^f{rf+4h{r)^^^ + 

+a(r)^ [256 + 48r2-3/ + 45r'^/(r)'^ + 2r2 (^-40 + 89r2) /z(r)^ - 115r4/z(r)'^+ 
+ 2rVW^ (^8 + 5 - 35 r2^(rf)] - 
-2r^a(r) [53r^/(r)^+4 (^-l+/z(r)^) (^-16 - 3r2 +7r^/i(r)^) -/W^ ^32+ 17r^ +47r^/i(r)^)] + 
[l6 + 3 + 3 /(r)^ - (32 + 19 r^) /i(r)^ + (^16 + 23 r^) /z(r)^ + 33 h{rf-\- 

+ /(r)4 (52-r2 + 39r2/i(r)2) +/(r)2 (^_24-5r2+ (-72 + 22r2) h{rf + 69 h{r)^)'\ } . 

tr now means that we have summed [c^jss and [c^]66- 

cf(r) = l-f(rf-h{rf 

c^{r) = ^ [3r2 + 2rVW^-(5r2+4a(r)2)/z(r)2 + 2r2/z(r)V 

+/(r)2 (^-4 - 5 + 8 a(r) - 4 oc(r)^ + 4 /z(r)^) ] 
cf{r) = ^|-10r4 + 4rVW^+ [23r4-8r2a(r) + 16 (l + r^) a(r)V 32 a(r)^ + 16 a^"^] /z(r)^+ 
[-17r2 + 8a(r)-16a(r)^] ^^"^ +4r4/z(r)^ + 

+r2/(r)'^ [-24 - 17 + 40 a(r)- 16 a(r)^ + 12 r2/i(rfj + 

+/(r)^ [-32a(r)^ + 16 a(r)'^ + 8 r^a(r) (^-5 + 6/i(r)^) + 16a(r)^ (^1 - 2r^/i(r)^) + 
+ (^24 + 23r2-2 (^10+ Hr^) /i(r)^ + 12r2^(r)'^)] | . 



3.6.2 One-loop SSTS mass shifts 

In this sub Section we offer some Tables with Mathematica calculations of the one-loop cor- 
rections to the masses of topological solitons in the semi-local Abelian Higgs model. We denote 
the Seeley coefficients, calculated by means of numerical integration of the circle symmetric Seeley 
densities evaluated at the numerical solutions for the field profiles, in the form: 

poo poo poo 

trci = 27t drrtrciir) , trc^ = 2n drrixc^^r) , c^ = 2n drrc^{r) 
Jo Jo Jo 

First, we focus on the embedded ANO vortices, the ho topological solitons. Second, we give results 
for several values of ho up to a value of ho close to /iq = 1, which corresponds to the CP^ -lumps. 
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1. One-loop vortex mass shifts in units of hm: 

Because the ANO vortex solutions were generated numerically, integration over the whole 
plane of the Seeley densities can also only be performed numerically. Therefore, we are forced 
to put a cut-off in the area and replace the infinite plane by a discus of radius R, which in the 
calculations shown below was chosen to be = 10.000 





ho = 0.0 


n 


trc' 


trc^ 




1 


-41.4469 


-91.8429 


12.599 


2 


30.3736 


0.96286 


2.61518 


3 


12.9447 


-0.0592415 


0.32005 


4 


4.22603 


0.001512548 


0.0230445 


5 


1.05059 


0.000758663 


0.0013023 


6 


0.20900 


-0.00023912 


0.0000698185 









/ = 1 


2 


-1.61536 


3 


-1.66862 


4 


-1.67809 


5 


-1.67966 


6 


-1.67989 



Table 7: Seeley coefficients for the (left) and Quantum Mass Correction (right) to the soliton in the semi- 
local Abelian Higgs model with /iq = 0.0. 



In the Table on the left, the first six Seeley coefficients are given. The next coefficients are 
very small and one expects that the approximation to the exact value of the one-loop mass shift 
is quite good, as shown in the Table on the right: the mass shift obtained by counting five or six 
coefficients agrees up to the third decimal figure. In fact, the thresholds in the wells are < 1 . 1 is the 
expectation value of the scalar field at the vacuum, and we find a similar situation to sine-Gordon 
kinks (see the first Part), where keeping six coefficients provides a fairly good approximation. 

From the same Table one can read the one-loop mass shifts of ANO vortices with / = 1 in the 
Abelian Higgs model: just take trc^(^) equal to zero and subtract two zero modes, not 4. The ratio 
is: 

1.67989 



^SLAHM 



1.09449 



= 1.53486 



A similar proportion exists between the ratios of kink mass shifts in the A (0)2 model and the BNRT 
model, both treated in the first Part: 



AM|™ 
AM^^M 



0.693943 
0.471113 



1.47299, a = 0.99 



AMfNRTM 
AM^^M 



0.698445 
0.471113 



1.48254, C7 = 1.01 . 



2. SSTS one-loop mass shifts 

Whereas the one-loop mass shift of ANO embedded vortices is always negative and varies 
extremely slowly as the area increases towards more negative values, one-loop mass shifts of gen- 
uine semilocal topological solitons with \ho\ > become less negative, and even positive, for larger 
areas, as is shown in the following Table. 

The classical degeneracy in energy between semi-local topological defects seems to be broken 
by one-loop fluctuations, the ANO embedded vortices becoming the ground states in the topolog- 
ical sector of one quantum of magnetic flux. It is remarkable how strong this effect becomes for 
topological solitons close to CP^-lumps, |/io| — 1. 
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R 


AMv{No = 6,R) 


AMv{No = 6,R) 


AMv{No = 6,R) 


AMv{No = 6,R) 


AMv{No = 6,R) 




ho 0.0 


ho = 0.1 


ho 0.3 


ho 0.6 


ho 0.9 


10^ 


-1.67955 


-1.61672 


-1.05000 


2.10142 


24.6066 


10^ 


-1.67971 


-1.58311 


-0.626167 


4.5485 


42.7747 




-1.67989 


-1.55133 


-0.252586 


6.41655 


60.9433 


10^ 


-1.68005 


-1.51957 


0.12086 


8.5741 


79.1116 


10^ 


-1.68026 


-1.48779 


0.49433 


10.7203 


97.2798 



Table 8: One-loop mass shifts for semi-local topological solitons: Five values of ho, five values of R, and 
fixed No = 6. 



3.6.3 Infrared divergences: quantum fate of semi-local topological solitons 

The origin of the degeneracy breaking is the slow decay (non-exponential) to their vacuum 
values of genuine semilocal topological solitons as compared with ANO vortices. Plugging the 
asymptotic form of the circle synmietric topological soliton solutions in the Seeley densities, we 
find the following behavior at infinity in terms of the parameter \h^\ (which sets the long r behavior 
of the solutions): 

27rrtrc{{r) _ p) + 4j(i2|/^i p _ |4) + ^(_L) 

2;rrtrcf(r) -^(i + + 4j(4|/,i|2 + |/,i|4) + ^(_L) 

Rtt 1 
2nrcf{r) _|/ji|2 + ^(^) 

27crtrc^{r) ^\h'\^ +^(-l+4\h'\^ -\h'f) + 
2Kr4{r) ^|/ji|2 + ^(768|/ji|2-80|/jir) + ^(l) 

2Kri,c%) -?^|/,i|2 + |!^(_4|/ji|2 + |,ji|4) + ^(^) 
2nrc-[r)'- '-^\h^? + ^{^) 
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2nn.cOir) ^^^^ + + 

The key observation is the appearance of infrared logarithmic divergences in the Seeley coefficients 
trc^(7^) and iYc^{K) for all n. The ghost coefficients c^{K^), however, are infrared convergent. The 
combination of the signs that we have seen in the previous sub-sections with the long r behavior 
shows that one-loop mass shifts of semi-local topological solitons tend to +00 in the infinite area 
limit. Semi-local topological defects grow infinitely massive due to the infrared effects of one-loop 
fluctuations. This phenomenon seems to be amazingly close to the non-existence of Goldstone 
bosons in (l+l)-dimensions. 

There is a very important exception: for ANO vortices, = and only the first-order co- 
efficients are infrared divergent. However, the contribution of these coefficients is totally canceled 
by mass renormalization counter-terms. Our results suggest that only the ANO vortices between 
all the semi-local topological solitons survive one-loop quantum fluctuations. It would be very 
interesting to try a more analytic approach to this problem in order to fully elucidate this delicate 
issue. 

3.6.4 Circle symmetric self-dual Abrikosov-Nielsen-Olesen vortices: one-loop mass shifts up 
to/ = 4 

Finally, we consider the problem of computing the one-loop mass shifts for superimposed 
ANO vortices at the /c = 1 limit up to four quanta -Stt- of magnetic flux in the Abelian Higgs 
model (no fluctuations in the a = 5,6 directions and 21 zero modes). The Figure 13 shows the 




Figure 13: Plots of the field profiles a(r) (a) and /(r) (b), the magnetic field B{r) (c), and the energy density 
e(r) for self-dual vortices with / = 1 (solid line), / = 2 (broken line), / = 3 (broken-dotted line) and / = 4 
(dotted line). 

field profiles, the magnetic field, and the energy density for the numerically generated solutions for 
/ = 1, / = 2, / = 3, and / = 4 in [29]. The Figure 14, however, encompass the 3D plots of these 
solutions. 

The Seeley coefficients for the different values of / as well as the one-loop mass shifts are 
displayed in the Tables below ^ 
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r 



Figure 14: 3D graphics of the energy density for / = 1, / = 2, / = 3 and / = 4 self-dual symmetric ANO 
vortices. 





l = \ 


1 = 


2 




/ 


= 3 


1 = 


4 


n 


ircn{K) 




trcn{K) 




n 


trc„(7i:) 




tVCn{K) 




2 


30.36316 


2.60773 


61.06679 


6.81760 


2 


90.20440 


11.51035 


118.67540 


16.46895 


3 


12.94926 


0.31851 


25.61572 


1.34209 


3 


36.68235 


2.60898 


46.01141 


4.00762 


4 


4.22814 


0.022887 


8.21053 


0.20481 


4 


11.69979 


0.46721 


14.64761 


0.77193 


5 


1.05116 


0.0011928 


2.02107 


0.023714 


5 


2.86756 


0.067279 


3.58906 


0.11747 


6 


0.20094 


0.00008803 


0.40233 


0.002212 


6 


0.566227 


0.0079269 


0.667202 


0.01620 



Table 9: Seeley Coefficients for / = 1, / = 2, / = 3 and / = 4 self-dual symmetric ANO vortices. 



No 


AMy(iVo) 


AMv{No) 


AMv{No) 


AMy(iVo) 




/ = 1 




/ = 3 


/ = 4 


2 


-1.02951 


-2.03787 


-3.01187 


-3.97025 


3 


-1.08323 


-2.14111 


-3.15680 


-4.14891 


4 


-1.09270 


-2.15913 


-3.18208 


-4.18014 


5 


-1.09427 


-2.16212 


-3.18628 


-4.18534 


6 


-1.09449 


-2.16257 


-3.18690 


-4.18606 



Table 10: One-loop mass shift for / = 1, / = 2, / = 3 and / = 4 self-dual symmetric ANO vortices. 



and the last Table provides the one-loop mass shifts of circle symmetric self-dual ANO vortices up 
to four quanta of magnetic flux taking into account A^o = 6 Seeley coefficients. 

^The Seeley densities corresponding to self-dual ANO vortices superposed at the origin up to four quanta of mag- 
netic flux are given in [29]. 
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/ 




AM 


1 


-1.09449 


-1 


2 


-2.16257 


-2 


3 


-3.18690 


-3 
-4 


4 


-4.18606 





It is remarkable that the one-loop mass shift seems to be linear in /. We have found that the one- 
loop correction to the mass of a circle symmetric vortex of magnetic flux 8 tt is extremely close 
to four times the one-loop mass shift of a similar vortex of magnetic flux 2n. This result strongly 
supports our method: after extremely sophisticated calculations we end with the natural answer. 

4. Prospects on the future of the subject 

We finish by offering a summary of other possible approaches to this subject as well as pointing 
to other playgrounds where similar methods may work and should be applied. 

1. In [49], and [50] Bias and Carrion explored a generalized sine-Gordon model containing 
as many scalar fields as the rank of the §L(A/^,C) Lie group. The potential energy density 
is determined from the simple roots in such a way that the integrability of the sine-Gordon 
model is enjoyed by the generalized model. Clearly, there is a correspondence between this 
model and the ordinary sine-Gordon model that is analogous to the correspondence between 
the deformed linear 0(A/^)-sigma model and the model. There are several types of 
kinks in this system, which have also been dealt with by the authors in a quantum setting. 
Nevertheless, it seems to us that the generalized sine-Gordon models offer another excellent 
arena for applying and improving the quantization method described in this report. 

2. In our Lectures we have addressed a massive non-linear S^-sigma model having a rich kink 
manifold, see [33]. In fact, we also computed the one-loop correction to the classical mass 
of the topological kinks using the Cahill-Comtet-Glauber formula in [34]. We look forward 
to attacking this problem by means of the heat kernel/zeta function approach, not only in this 
system but also in other similar non-linear sigma models where we vary the potential energy 
density; e.g., adding quartic terms, and/or we change the target space, e.g., to 

3. In a very interesting paper [51] Vafa et al. analyzed ^ = 2 supersymmetric Landau- 
Ginzburg models that are integrable deformations of the ^ = 2 supersymmetric minimal 
series. The BPS states of these integrable (1 + 1) -dimensional field theories were identified 
by these authors as "holomorphic" kinks with very noticeable properties. Almost ten years 
later several of us went through the same set of kinks in [52]. We focused only upon the 
bosonic sector of the model and chose a real analytic point of view. We believe that the 
techniques developed in these Lectures provide a procedure for computing the one-loop kink 
mass shifts due only to bosonic fluctuations, is the threshold of the wells and hence we 
expect to need at least A/q = 20 coefficients in order to reach a good approximation. In the 
o/K = 2 supersymmetric version of this model, the one-loop correction including bosonic and 
fermionic fluctuations is zero. Therefore, there is no need to compute the effect of fermionic 
fluctuations: bosons and fermions created in the kink background cancel each other exactly. 
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4. Another strategy may be fruitful. The heat trace, or partition function, which is basic in 
our approach, admits another conceptual understanding as a path integral over closed world 
lines: 

1 ^ dZa dZa 



Se[z\,--- ,zn] ^ I dT 
Jo 

Z{P,l) = Tr^2^-^^ = n^^i dxa r^'^"'''^z(T)^-^^^[^i'-'^^l . (4.1) 



Here, T = it is imaginary time, 5^ is the Euclidean action for a particle moving in a cube 
\^ with varying positions z^?(t). The path integral is over all the closed paths, and ordinary 
integration over all the base points of the loops is necessary. It is well known that the Feyn- 
man treatment of these path integrals in the high-temperature regime [53] reproduces the heat 
kernel expansion. What seems more promising is the numerical computation of these world 
line path integrals using Montecarlo methods, see e.g [54]. Applied to the ANO vortices and 
semilocal strings this should provide reliable results to be contrasted against our calculations. 



5. The Lagrangian density 



/ 11 1 cosh (7Z 



^\x^) (4.2) 



describes the dynamics of a vara avis scalar quantum field theory. The mesons do not move 
freely even though they do not interact because they are constrained by a background. If 
(7 = 0, the background in (4.2) is formed by two parallel plates located at a distance / from 
each other in two planes orthogonal to the z-axis. The effect of the plates is mimicked by two 
5 -function potentials of strength A±. When / = oo and O" = 1, the background corresponds 
to a kink living in the z-axis, or to a solitonic/thick domain wall orthogonal to the z-axis. 
Other choices of the couplings provide more complex backgrounds built from these two 
basic backgrounds. 

In Reference [55] Milton wrote the Green's function for the two-plate setup (using Dirichlet 
boundary conditions). This leads to the energy momentum tensor encoding the Casimir 
energy -essentially the Too component- and the Casimir force -essentially the T^z component-. 
It is tempting to perform the same calculation for the kink background. This should provide 
not only information about the one-loop kink mass shift but should also shed light on the 
qualitative nature of the forces exerted by the scalar fluctuations on the kink profiles. 

6. Over the last year, an interesting paper by Baacke and Kevlishvili was published [56] in 
which the one-loop shifts to the classical masses of Nielsen-Olesen vortices were obtained 
by using Green's function methods. The remarkable fact is that the authors gave the quantum 
corrections with no restriction in the ratio of Higgs and particle masses. 

It seems that the time is ready to tackle the problem of the quantum corrections of Z- 
electroweak strings. These topological defects are embedded NO vortices in the neutral 
massive vector field of electroweak theory, and the calculation may will be of experi- 
mental interest (even though some imaginary contribution to the energy will arise, because 
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electroweak strings are unstable). In fact, there is an interesting paper on this subject, see 
[57], at the non-physical value of the weak angle equal to zero, and even a status report [58]. 

7. In a longer perspective, one might think of studying the quantum fluctuations of BPS mag- 
netic monopoles. The work in [59] suggests that the bosonic sector of ^ = 2 supersymmet- 
ric Yang-Mills is the right model to look into this problem. Even though van Nieuwenhuizen 
et al. succeeded in computing the one-loop mass shift to ^ = 2 SUSY monopoles, the diffi- 
culties in a purely bosonic framework seem insurmountable. First, the second-order operator 
governing the fluctuations is (for the SC/(2) group) a 21 x 21 matrix-Schrodinger operator 
in three dimensions. Second, the theory is renormalizable, not super-renormalizable. The 
coupling constant also receives one-loop divergent contributions that must be canceled by 
the secod-order Seeley coefficient. So who is afraid of this big bad wolf? 
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